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PREFACE 

WHILE  lecturing  in  1913  on  the  calculation  of  self  and  mutual 
induction,  the  writer  noticed  that  the  integral  in  the  case  of  co- 
axial circles  was  of  such  a  form  as  to  make  the  application  of  the 
arithmetico-geometrical  means  obvious.  An  extremely  elegant  formula 
(Proc.  Roy.  Soc.  A,  Vol.  c.  1921,  p.  63.  Appendix,  this  book,  No.  35) 
when  put  to  numerical  test  justified  the  method  as  one  well  adapted  to 
computation.  The  calculation  of  more  complex  cases  of  self  and  mutual 
induction  led  to  several  interesting  formulae  for  the  complete  third 
elliptic  integral.  With  a  view  to  making  a  complete  study  of  the  group 
of  recurrence  formulae  associated  with  the  A.G.M.  scales,  further  work 
was  postponed  until  a  vacation  in  1921,  when  many  of  the  results  in- 
cluded in  this  book  were  obtained,  and  were  thought  by  the  writer  to 
be  new.  On  later  consulting  the  Collected  Works  of  Legendre,  Gauss 
and  Jacobi,  the  writer  found  that  he  had  rediscovered  many  known 
formulae,  but  that  quite  a  number,  likely  to  be  of  service  in  computation, 
appeared  to  be  hitherto  unknown.  In  the  summer  of  1922  the  writer 
decided,  as  a  vacation  task,  to  include  in  a  single  monograph  the  entire 
theory  of  elliptic  functions  associated  with  the  A.G.M.  scales,  thus  com- 
pleting, and  to  some  extent  adding  to,  the  work  of  Gauss  left  unpublished 
after  his  death,  and  placing  before  mathematicians  in  accessible  form 
a  mode  of  approach  to  Elliptic  Function  Theory  directly  related  to  the 
art  of  machine-computation.  It  is  interesting  to  note  in  this  connection 
that  this  subject  was  to  have  formed  the  content  of  the  third  volume 
of  Halphen's  Functions  Elliptiques,  unfortunately  left  incomplete  at  the 
time  of  his  death. 

While  the  present  monograph  contains  many  new  formulae  and 
methods  of  computation,  no  claim  is  made  as  to  novelty  in  fundamental 
analytical  treatment.  References  to  more  important  books  and  memoirs 
consulted  by  the  author  are  given  in  foot-notes.  A  complete  bibliography 
is  not  given  since  the  reader  will  find  exhaustive  references  in  the  Royal 
Society  Index,  1800-1900,  Vol.  I,  Mathematics,  and  after  1900  in  the 
mathematical  volumes  of  the  International  Catalogue  of  Scientific 
Literature. 
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ON  THE  DIRECT  NUMERICAL  CALCULATION 
OF   ELLIPTIC    FUNCTIONS   AND   INTEGRALS 

SECTION  I 
INTRODUCTION 

ALTHOUGH  the  numerical  evaluation  of  elliptic  functions  and  integrals 
is  dealt  with  in  many  of  the  standard  treatises,  it  cannot  be  said  that 
the  use  of  tables  is  entirely  satisfactory  in  dealing  with  computations 
required  in  several  branches  of  physics  and  astronomy.  Many  of  the 
integrals  appearing  in  these  problems  can  only  be  expressed  as  some- 
what complicated  expressions  involving  the  complete  elliptic  integrals 
of  the  first  and  second  kinds,  and  the  corresponding  functions  of  com- 
plementary modulus.  These  formulae  give  rise  in  many  cases  to 
differences  of  nearly  equal  quantities :  the  numerical  calculation  of  such 
expressions  to  a  sufficient  number  of  significant  figures  is  thus  often  a 
matter  of  some  difficulty,  requiring  tedious  and  sometimes  uncertain 
interpolations. 

While  the  numerical  calculation  of  the  Theta-functions  is  extremely 
rapid  by  the  use  of  the  highly  convergent  series  involving  the  name  q, 
these  are  not  usually  the  functions  which  make  their  appearance  directly 
in  physical  or  astronomical  problems.  The  necessary  transformations 
required  to  express  the  results  in  terms  of  the  (/-series  require  in  many 
cases  somewhat  complicated  analysis,  while  numerical  computation  by 
this  method  necessitates  the  exact  calculation  of  a  large  number  of 
auxiliary  quantities*. 

The  forthcoming  tables  in  course  of  construction  under  the  auspices 
of  the  British  Association  t  are  so  arranged  as  to  allow  a  large  number 
of  functions  to  be  evaluated,  and  will  very  greatly  aid  numerical  work 
in  many  practical  problems  involving  the  use  of  elliptic  integrals  and 
functions.  While  sufficient  for  most  purposes,  the  use  of  these  tables 
will  still  involve  considerable  interpolation  when  it  is  required  to  take 
out  values  corresponding  to  arguments  and  moduli  expressed  by  six  or 

*  On  this  point  see  Eosa,  E.  H.,  and  Grover,  F.  N.,  Bulletin  169,  Bureau  of  Stan- 
dards, Washington  (1913),  pp.  67  and  73;  Nagaoka,  H.,  Tok.  Coll.  Sc.  J.  27  (1909), 
No.  6,  and  Jour.  Tok.  Math.  Phys.  Soc.  4  (1908),  p.  284;  6  (1911),  p.  10;  Nagaoka  and 
Sakurai,  Scientific  Papers  of  the  Institute  of  Physical  and  Chemical  Research,  Vol.  n. 
Dec.  1922,  pp.  1-67;  Olshausen,  G.  K.,  Phys.  Rev.  51  (Dec.  1910),  pp.  617-636. 

t  See  report  (by  Sir  A.  G.  GreenhiU),  Brit.  Ass.  Report,  1912  (Dundee),  pp. 
39-55. 
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seven  significant  figures,  as  is  now  necessary  in  several  modern  applica- 
tions of  elliptic  functions. 

For  several  reasons  it  is  advantageous  to  have  available  a  direct 
method  of  computation,  independent  of  auxiliary  tables,  based  on  the 
use  of  the  modern  high-speed  calculating  machine.  The  writer  has 
therefore  thought  it  desirable  to  set  out  in  detail  the  method  of 
calculation  based  on  Landen's  quadric  transformation  and  the  use  of 
the  scale  of  arithmetico-geometrical  means,  recapitulating  briefly  the 
original  developments  of  Lagrange,  Legendre,  Gauss,  Jacobi  and  others, 
extending  them  in  several  directions  and  keeping  in  mind  the  ultimate 
use  of  the  calculating  machine  for  numerical  work. 


SECTION  II 

HISTORICAL  NOTE  ON  LANDEN'S  TRANSFORMATION  AND 
THE  VARIOUS  SCALES  OF  MODULI  AND  AMPLITUDES 

The  importance  of  realizing  rapid  and  accurate  methods  of  calculating 
the  elliptic  integrals  now  denoted  by 

u  =  F(<l>,  ti)=t  <fy/A U,  k}  and  EU,  k)  =  /*Afo,  k)  <fy  ...(1), 
h  Jo 

where  A  (<£,  K)  =  V(l  -  #  sin2  <£), 

was  first  remarked  by  Euler*  (1766),  although  it  was  not  until  several 
years  later  that  Landent  (1775)  discovered  in  geometrical  form  the 
transformation  which  forms  the  basis  of  the  existing  methods  of  the 
numerical  calculation  of  the  elliptic  integrals. 

A  method  of  successive  transformations  for  the  ultimate  reduction  of 
the  algebraic  forms  of  these  integrals  to  elementary  integrals  was  pub- 
lished by  Lagrange  J  in  1784-85.  This  memoir  contains  an  exposition  of 
the  scales  of  arithmetico-geometrical  means,  together  with  two  types  of 
algebraic  transformations  corresponding  to  the  increasing  and  diminishing 
amplitudes  of  Landen's  trigonometrical  forms:  it  is  shown  that  the 
given  integral  must  lie  between  two  limits  which  may  be  made  to 
approach  each  other  as  closely  as  may  be  desired :  the  limits  thus  ob- 

*  Euler,  Novi  Comm.  Acad.  Sc.  Petrop.  vol.  x.  1766. 

t  Landen,  Phil.  Trans.  Roy.  Soc.  vol.  LXV.  p.  283,  1775;  Mathematical  Memoirs, 
London,  1780. 

J  Lagrange,  "Sur  une  nouvelle  methode  de  calcul  integral  pour  les  diff^rentiellea 
affect^es  d'un  radical  carr£  sous  lequel  la  variable  ne  passe  pas  le  quatrieme  degre," 
Mem.  de  VAcad.  roy.  des  Sc.  de  Turin,  t.  n.  1784-5;  (Euvres  (Gauthier- Villars,  Paris, 
1868),  t.  n.  pp.  253-312. 
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tained  for  the  integrals  of  the  first  and  second  kinds  are  practically 
identical  with  the  limits  and  series  obtained  about  the  same  time  by 
Legendre. 

Landen's  transformation  was  applied  by  Legendre*  to  the  numerical 
calculation  of  the  elliptic  integrals,  and  furnished  the  method  of  com- 
putation by  means  of  which  the  latter's  well-known  tables  were  con- 
structed. The  scales  of  moduli  and  amplitudes  employed  by  Legendre 
are  briefly  described  below. 

If  we  form  successive  moduli  and  amplitudes  according  to  the 
recurrence  formulae 

#,.,.!  =  (1  -*.')/(!+*«')  «»d  tan(«k+1-<k)  =  *»'tan«k  -.-(2), 
Landen's  transformation  in  trigonometrical  form  leads  to  the  result 


As  n  increases,  the  successive  moduli  kn  converge  to  zero,  and  the 
complementary  moduli  &„'  to  unity,  while  the  amplitudes  increase  so 
that  <£n/2n  tends  rapidly  to  a  definite  limit  :  \iltimately  we  may  write 
F(<j>n,  £„)  =  <£»  and  thus  obtain  the  result 

F(4*,  *0)  =  (a*/*-)  6,/2»,  where  K=  i»  (1  +  *0  (1  +  *0  ...  (1  +  *«) 

.........  (4). 

The  result  of  calculating  successive  moduli  and  amplitudes  in  reverse 
order  is  equivalent  to  calculating  in  forward  order  these  quantities 
according  to  the  recurrence  formulae 

kn+1  =  2jknl(l  +  k^)  and  sin  (2^,,+!-^)  =  *,,  sin  ^«   ...(5). 
This  transformation  then  gives 

-F(*.«,*.«)=Ki+*0^(*.,*0  ...............  (6). 

As  n  increases  the  successive  moduli  kH+i  converge  to  unity  and  the 
amplitudes  to  a  limit  <^B  :  we  thus  obtain,  when  n  is  sufficiently  large, 


sec  <M<A»  =  log  tan  (%*  +  Jfc.)    ......  (7), 

giving  finally 

J?Oo,  *•)  =  V(*itt  -  40A/V  logtan  frr  +  J^,)     ...(8). 
Legendre  also  developed  for  purposes  of  numerical  calculation  extremely 
convergent  series  for  E($,  k)  which  he  employed  in  the  tabulation  of 
this  function. 

*  Legendre,  "  Memoire  sur  les  integrations  par  arcs  d'ellipse"  and  "Second 
memoire,  etc.,"  Mem.  de  VAcad.  des  Sciences  de  Paris,  ann.  1786  (Paris,  1788), 
pp.  616-643  and  644-683;  Traite  des  Fonctions  Elliptiques,  Paris,  1825,  t.  L 
p.  79  et  seq. 

1—2 
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The  method  discovered  independently  by  Gauss  and  called  by  him 
the  "algorithm  of  the  arithmetico-geometrical  mean"  originated  in 
connection  with  the  evaluation  of  complete  elliptic  integrals  arising 
from  a  problem  in  attractions  required  in  planetary  theory*:  in  a  note  he 
mentions  that  the  results  were  obtained  by  him  as  part  of  a  more  com- 
prehensive theory,  independently  of  the  results  of  Lagrange  and  Legendre 
with  which  they  are  closely  connected.  Gauss  employed  a  trigono- 
metrical transformation  along  the  following  lines. 

Commencing  with  two  numbers  (a0,  60),  successive  numbers  («„,  bn) 
are  calculated  from  the  recurrence  formulae, 

«n+l  =  K«n  +  bn),     &n  +  l  =  V(«A)    ...............  (9). 

Writing  A^  =  >/(«w2  cos2  <f>n  +  bn2  sin2  <£n),  Gauss  employs  the  recurrence 
formula 

an  tan  <f>n+1  =  An+1  tan  <£n, 

or  tan  <£n+1  =  tan  <f>n  ,J(an+lla^  .  J(bn  +  An)/  J(an  +  An)  .  .  .(W). 

From  the  relation 

tan<£n  =  (AaA2...  An...)/(>0aa  ...  «„...).  tan  <£0 
we  have  ultimately,  since  the  <£'s  tend  to  a  limit  <£„, 


=  <£„/«„, 
o 
and  hence 

«»  I  °  ^o/A0  =  tan-1  [(A!  A2...An...  )/(«„«!...«„...).  tan  <#>„]  ...(11). 
./o 

Extremely  convergent  series  for  the  elliptic  integral  of  the  second 
kind  were  also  obtained  by  Gauss.  The  application  of  this  method  to 
the  numerical  calculation  of  the  elliptic  integrals  was  considered  by 
Jacobi  some  years  later  t. 

In  the  course  of  researches  carried  out  between  1797  and  1808,  but 
not  published  until  after  his  death  considerably  later  J  ,  Gauss  developed 
the  theory  of  elliptic  functions  to  a  remarkable  extent  and  derived  a 
large  number  of  formulae  which,  at  a  later  date,  were  obtained  inde- 
pendently by  Abel  and  Jacobi  and  which  are  now  associated  with  the 

*  Gauss,  "  Determinatio  Attractionis  quam  in  Punctum  quodvis  Positionis  datae 
exerceret  Planeta,  etc.,"  Comm.  Gott.  Soc.  Eeg,  Scient.  iv.  1818;  Ges.  Werke 
(Gottingen,  1866),  Bd.  in.  pp.  331-355.  The  note  mentioned  above  is  dated  Feb.  9, 
1818  ;  Ges.  Werke,  Bd.  in.  pp.  357-360. 

t  Jacobi,  Fundamenta  Nova,  1829,  sections  38  and  52  ;  Ges.  Werke  (Berlin,  1881), 
Bd.  i.  pp.  154  and  203.  Also  "  Numerische  Berechnung  der  Elliptischen  Functionen," 
Crelle,  Bd.  xxvi.  pp.  93-114;  Ges.  Werke  (1861),  Bd.  i.  pp.  345-368. 

+  Gauss,  Ges.  Werke  (1866),  Bd.  in.  pp.  361-403  and  433-469. 
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theory  of  the  Theta-funetions.  In  one  of  these  posthumous  papers,  Gauss* 
employs  the  recurrence  formula  tan2<£n  =  >/(6»-n/am-i)tan<^B+i,  which 
is  much  more  convenient  for  purposes  of  numerical  calculation  than  (10). 
The  last  mentioned  method  is  considerably  simpler  than  that  of 
Legendre's  scales  of  moduli,  especially  when  a  calculating  machine  is 
available,  although  in  theory  these  two  methods  of  computation  are 
practically  identical. 

Although  a  few  formulae  based  on  the  scales  of  arithmetico-geometrical 
means  were  briefly  considered  by  Jacobit,  no  attempt  seems  to  have  been 
made  by  him  to  apply  the  results  to  actual  computation. 

As  the  method  of  evaluating  the  elliptic  functions  by  the  use  of  the 
A.G.M.  scales  is  only  briefly  dealt  with  in  existing  text-books  J,  and  does 
not  seem  to  have  been  developed  to  the  extent  it  deserves  in  various 
memoirs  which  have  appeared  on  the  subject  in  recent  years,  a  tolerably 
full  account  of  the  subject  is  given  in  the  present  book,  together  with 
a  brief  outline  of  the  main  analytical  steps.  Several  new  formulae  have 
been  obtained  by  the  writer  which  facilitate  the  application  of  this 
-method  of  numerical  computation  and  considerably  extend  its  scope. 

From  the  point  of  view  of  pure  analysis,  the  modern  "Weierstrassian 
notation  presents  many  advantages  in  regard  to  elegance  and  symmetry, 
a  superiority  which  no  longer  holds,  however,  when  it  is  required  to 
proceed  to  numerical  evaluation.  As  the  main  point  of  the  present 
paper  is  concerned  with  improvements  in  methods  of  computation,  the 
older  Jacobian  notation  has  been  adhered  to  as  far  as  possible,  especially 
as  it  expresses  most  readily  the  elliptic  functions  and  integrals  which 
arise  from  integrations  involving  circular  and  hyperbolic  functions.  The 
results  are  set  out  in  such  a  manner  as  to  provide  a  fairly  complete 
compendium  of  formulae  intended  to  be  of  service  to  the  computer  in 
reducing  to  numerical  results  such  expressions  involving  elliptic  functions 
and  integrals  as  may  arise  in  physical  and  astronomical  problems. 

*  Gauss,  loc.  cit.  p.  388. 

t  Jacobi,  Ges.  Werke  (1881),  p.  357. 

J  Cayley,  A.,  Elliptic  Functions,  2nd  Ed.  (G.  Bell  and  Sons,  London,  1895), 
pp.  326-338;  Enneper-Muller,  Elliptische  Funetionen  (2nd  ed.  1890),  pp.  361-364; 
Tannery-Molk,  Fonctions  Elliptiques,  t.  rv.  1902,  pp.  269-275.  According  to  a  foot- 
note in  Halphen's  Fonctions  Elliptiques,  i.  n.  1888,  p.  310,  the  subject  of  the 
arithmetico-geometrical  mean  was  to  have  received  special  attention  in  the  author's 
third  volume,  unfortunately  incomplete  at  the  time  of  his  death. 
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SECTION  III 

ON  THE  SCALE  OF  ARITHMETICO-GEOMETRICAL  MEANS 

In  forming  the  scale  of  arithmetico-geometrical  means,  we  start  with 
the  positive  numbers  (a0,  b0),  of  which  «0  is  the  greater,  and  form 
successively 


...(12). 

•  >         ' 

We  obviously  have  the  relations 


=       — 


cn  = 


The  a's  and  b's  tend  to  the  same  limit  denoted  by 

M(a0,bQ')  =  lim  an  ........................  (14) 

n-*-o> 

with  extraordinary  rapidity,  even  when  a0  and  b0  are  initially  numbers 
of  very  different  magnitudes*. 

We  notice  in  passing  that  M(a0,  60)  satisfies  a  homogeneity  relation 
of  the  form 

fM(a0,b0)  =  M((a0,eb0')  .....................  (15), 

e  being  any  number.  The  array  of  numbers  (12)  will  be  referred  to  as 
the  "  scale  of  arithmetico-geometrical  means  («„,  bo),"  or  more  briefly 
as  the  A.Q.M.  scale  (a0,  £><>)•  Also,  the  limit  M(a0,  b0~)  will  be  denoted 
by  an,  as  long  as  by  so  doing  no  ambiguity  is  involved. 

If  we  calculate  in  the  same  way  an  array  of  numbers  commencing 
with  {a0'  =  a0,  b0'  =  c0=  */(a<?-b02}},  we  derive  what  will  be  called  the 
"complementary  A.G.M.  scale  («„',  £<>')»"  the  symbols  being  in  this  case 
denoted  by  accented  letters,  and  the  limit  M(a0',  &„')  by  an'. 

*  The  usual  arithmetical  process  of  extracting  square  roots  is  easily  adapted  to 
most  calculating  machines  of  modern  type.  An  account  of  these  is  given  in 
Modern  Instruments  and  Methods  of  Calculation,  edited  by  E.  M.  Horseburgh 
(London,  G.  Bell  and  Sons,  1914)  ;  also  by  d'Ocagne,  Le  Calcul  Simplifig  (Gauthier- 
Villars,  Paris,  1905).  In  carrying  out  the  process  the  work  is  considerably  shortened 
by  keeping  in  mind  the  rule  that  if  the  first  p  digits  out  of  the  number  n 
required  in  the  square  root  have  been  obtained  by  the  usual  process,  the  next 
p  -  1  digits  can  be  obtained  by  division  only,  with  a  possible  error  of  1  in  the  last  digit 
(Chrystal's  Algebra,  Part  i.  5th  ed.  1904,  p.  210).  To  illustrate  the  rapid  convergence 
of  the  A.G.M.  scales,  Gauss  works  out  an  example  starting  with  a0  =  l  and  60  =  0'2  : 
<i5  and  ftfi  are  represented  by  the  same  number  to  fifteen  significant  figures. 
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It  is  obvious  that  the  operations  of  the  A.G.M.  scale  may  be  carried 
out  backwards  :  if  the  entries  are  denoted  by  negative  suffixes,  they 
are  connected  with  those  of  the  complementary  scale  by  the  relations 
a_B  =  2"aB',        b-n  =  2»cn',        c_B  =  2"6B' (16). 

These  relations  evidently  hold  good  when  accented  and  unaccented 
symbols  are  interchanged.  As  a  result  of  these  relations  it  will  be 
seen  that  there  is  no  advantage  to  be  gained  from  considering  any  other 
than  the  A.G.M.  scale  of  positive  suffixes  (aB,  &„)  and  the  corresponding 
complementary  scale  (an',  £>„'). 


SECTION  IV 

LANDEN'S  SCALE  OF  INCREASING  AMPLITUDES 

(i)   Calculation  of  the  Elliptic  Integrals  of  the  First  and  Second 
kinds,  modulus  k. 

If  we  write  A,  =  J(a^  cos2<£B  +  bn2  sin2«£B),  the  recurrence  formula 

tan(^»+i-^)  =  (^K)tan^»   (17) 

may  be  written  in  either  of  the  forms 

sin  (2<£n  -  <£B+1)  =  (C».H/«»+I)  sin  4>n+l, 

cos  (2<f>n  -  <#>B+1)  =  Aw+1/aB+1   (18). 

From  these  we  derive 

AB+i  +  cB+1  cos  <£„+!  =  An     and     AB+1-cB+1  cos  </>„+!  =  a,  &n/AB  ...(19), 
which  give  on  differentiation, 

(ty0/A0  =  Jd^/A1  =  id^Aa=...=(i)\fy^All  = (20). 

From  the  relation 

A2B+i  +  cn+1  An+1  cos  <£n+1  =  \  AB2  +  \  anb% 
we  obtain,  making  use  of  (20), 

&nd<t>n - &n+ld<f>n+l  =  cn+1  cos <t>n+ld<f>n+1  - aAd<£B/AB  --.(21). 
Subtracting  from  each  side  the  identity 

a«2<ty»/A»  -  «2B+i^«+1/An+1  =  Jc.'dk/A,  -  an  bnd<j>nj*n . .  .(22) 
we  obtain 
(An  -  aB2/An)^B  -  (AB+1  -  a2n+1/AB+1)^B+1 

=  cB+1cos<^)n+1^n+1-icB2^n/An  ...(23). 

Writing  down  a  series  of  such  equations,  commencing  with  n  =  0, 
and  noting  that  as  n  increases,  (AB  -  a»2/AB)  tends  to  zero  in  the  limit, 
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we  have,  on  integrating, 

fbo  /"<fro 

I    A0o?<£0  —  a02  I     d<f>0/&0  =  d  sin  <^  +  c2  sin  <£2  +  ...  +  cn  sin  <bn  +  ... 

^o  yo 

/vfc. 

, (24). 


If  we  now  construct  the  A.G.M.  scale  (a0  =  l,  b0  =  k'},  and  calculate 
successive  amplitudes  commencing  with  <£0  =  <£,  we  have  A0  =  A  (<£,  k). 
As  n  increases,  the  a's  and  b's  tend  to  the  same  limit  an,  while  the 
angle  <f>n/2n  tends  to  a  finite  limit.  We  thus  obtain  from  (20)  and  (24), 

w  =  l^fc*)  =  (1/0.)  to/2-)  ..................  (25), 

E((f>,  k)-F(<t>,  k}  =  (<?!  sin  <£x  +  c^  sin  <£2  +  •  ••  +  cnsm<f>n  +  ...) 

-Hc02  +  2c12  +  4c22+...+2^2+...).F(<M)  ......  (26). 

If  we  commence  with  <£0  =  |TT,  we  have 

&  =  T,  <£2  =  27r,  <£3  =  47T  .  .  .  <£n/2n  =  Jir, 
so  that  K=\Tr\an,  and 

(^-JEr)/^=Kco2  +  2c12  +  4c22+...  +  2"cri2+...)  ......  (27).' 

The  series  (24)  and  (27)  converge  with  extreme  rapidity,  resembling 
in  this  respect  the  g-series  with  which  they  are  closely  connected*. 
Jacobi's  integral  Z  (<£,  K)  is  defined  by  the  relation 

Z(+,  k}  =  E(t,  K)-(EIIC)F(<}>,  *)  ............  (28) 

or,  if  the  angle  <f>  is  connected  with  the  argument  u  by  the  relation 

sin  <£  =  sn  (u,  k), 


Z(u,k}=     Anz(u,k')du-(E/K)u  ............  (29). 

^o 

From  (26)  and  (27)  we  derive  the  expansion 

Z(<j>,  K}  =  G!  sin  <#»!  +  c2  sin  ^>2  +  ...  +cTCsin<£n+  ......  (30). 

In  terms  of  Jacobi's  Theta-function  ®  (u,  k\ 

)  =  log  [0  («,  *)/©  (0,  *)]  .  .  .(31) 


o 
where  ®(0,X)  =  J(2VK/ir)   .....................  (32). 

*  These  relations  are  discussed  by  Legendre  (Traite  des  Fonctions  Elliptiques, 
t.  in.  Deuxieme  Supplement,  1828,  p.  Ill  et  seq.).  On  the  connection  between  the 
A.G.M.  scale  and  the  modern  transformation  theory  of  the  Theta-  functions,  see 
Tannery  and  Molk,  Fonctions  Elliptiques,  t.  iv.  1902,  note  2,  pp.  269-273.  A  develop- 
ment of  the  theory  of  the  Theta-f  unctions  based  on  Legendre's  scales  of  moduli  is 
outlined  by  Eichelot,  Correspondenz  mit  Herrn  Professor  Schroter  (Konigsberg, 
1868).  See  Examples  19-24  of  the  Appendix. 


f 

/• 
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From  (20)  we  have  d<f>/±(<t>,  yt)  =  (|)nc?0m/A,,  so  that  the  series  (30) 
enables  us  to  write 

r»  /•#!» 

Z  (<f>,  k")  <^>/A  (<f>,  k)  =  2  (cji/2*)  /     sin 
i  Jo 

On  integrating  and  making  use  of  (19)  and  (18),  we  obtain 

sin  fad<f>n/&n  =  (l/e,)  log  (a,-^.-!)  =  (l/cw)  log  sec  (2<£B_2  -  </>»-i), 
/o 

leading  to  the  expansions* 

f*  °° 

•M)  i 

or 

log  [®  («,  *)/0  (0,  *)]  =  |  log  (a0/A0)  +  i  log  sec  (2<£0  -  &) 

+ 1-  log  sec  (2<k  -  <fcj)  +  •  •  •  +  Q)**1  log  sec  (2<f>n_i  -  <£„)  4- (34). 

In  the  series  (34)  it  is  more  convenient  to  write 

for  purposes  of  numerical  calculation. 

'    (ii)  Calculation  of  the  Jacobian  Functions  sn  (u,  k),  cn  (u,  k),  dn  (u,  k) 
'    etc.  in  terms  of  the  argument  u. 

When  the  argument  u  is  given  and  the  modulus  k  is  known,  it  is 
only  necessary  to  compute  the  A.G.M.  scale  (a0=l,  b0  =  k'}  to  such  a 
value  of  n  that  cn  is  less  than  the  small  quantity  determining  the  order 
of  accuracy  of  the  calculations.  Since  u  =  F(<t>,  k),  equation  (25)  enables 
us  to  calculate  in  circular  measure  the  angle  <f>n  from  the  relation 

<f>n  =  2nanu   (35). 

From  the  recurrence  formula  (18), 

sin  (2<£n_i  —  <£„)  =  (Ca/#n)  sin  <£„ , 

.  we  are  enabled  to  calculate  successively  the  angles  <£„_!  ...<h,<J>i  and 
finally  <£0,  in  terms  of  which 
sn  (u,k)  =  sin  <l>0 ,   cn  (u,  k)  =  cos  <£0,   dn  (u,  k)  =  A0  =  cos  </>0/cos  (&  —  <£0) 

(36). 

Z(u,  k}  may  then  be  computed  from  the  series  (30),  and  ®  (u,  k)  from 
formulae  (32)  and  (34). 

*  Formulae  (30),  (33)  and  (34)  are  practically  the  same  as  those  given  by  Legendre 
(Traite  de»  Functions  Elliptiques,  t.  i.  section  90),  and  are  contained  implicitly  in 
formulae  given  in  Gauss'  posthumous  papers.  Formula  (34)  is  equivalent  to  that 
given  by  Jacobi  (Crelle,  Bd.  xxvi.  pp.  93-114;  Gesammelte  Werke  (1881),  Bd.  n. 
p.  357). 
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SECTION  V 

THE  HYPERBOLIC  SCALE  OP  INCREASING  AMPLITUDES 

(i)   Calculation  of  the  Elliptic  Integrals  of  the  First  and  Second 
kinds,  modulus  k'. 

Writing  i<f>n  for  <f>n  throughout  the  formulae  (17)  to  (24),  we  obtain, 

instead  of  (17), 

tanh  (<£n+1  -  <£„)  =  (bn/an)  tanh  <f>n (37), 

which  may  also  be  written 

sinh  (2<}>n  -  0n+1)  =  (cn+1/«n+i)  sinh  <£n+1 ,  cosh  (2<f>n  -  4>n+1)  =  ATO+1/an+1 

(38). 

If  we  denote  Are  =  <J(anz  cosh2  <f>n  -  bn  sinh2  <£„),  it  is  easily  proved,  as 

before,  that 

^0/A0  =  |^1/A1  =  i^2/A2=...  =  (|)«^n/An= (39) 

from  which  it  follows  that  as  n  increases,  we  have  in  the  limit, 

f  °<%/A0  =  (lK)(4>n/2n). 
Jo 
If  we  now  write 

sinh  <£0  =  tan  <f>    (40) 

it  is  easily  proved  that 

A0  =  A(<£,  £')/cos  <f>,     <tyo/A0  =  ^/A  (<£,  k'} (41), 

so  that 

u  =  F(4>)k'}  =  (\lan}(<t>nW)     (42). 

Formula  (24)  becomes, 

r<t>o  r<t>o 

AQ^O - «o2  I     ^0/Ao  =  (ci sinh *^i  +  c2 sinh <f>2+  ...  +cn smh <j>n  +  • . • ) 
Jo  Jo 

-Kco2  +  2c12  +  4c22+...+2*cn2+...)f  °^0/A0  ...(43). 

Jo 
It  is  easily  shown  that 

f  °A0^0=(  A(^>,^)^/cos2^ 
Jo  Jo 

=  I "dn2  (u,  k')  du/crf  (u,  k'}  =  I "dn2  (iu,  k}  du, 
Jo  Jo 

or  in  terms  of  the  ^-function  of  imaginary  argument, 

f    -Mfc  =  (\li}Z(iu,  k}  +  (EjK} u. 
Jo 

Substituting  in  (43),  and  making  use  of  (27),  we  derive  the  series, 
corresponding  to  that  of  (30), 

( 1  /i)  Z(iu,  £)  =  Cj  sinh  &  +  c2  sinh  <f>t  +  . . .  +  cn  sinh  <£„+...  ...  (44). 
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The  following  formula,  corresponding  to  (34),  easily  follows  by 
integration, 

log  [®  (tX  t)/&  (0,  *)]  =  |  log  (aj\)  +  i  (i)-+1  log  sech  (2^  -  *.) 

......  (45) 

in  which  we  may  conveniently  write  aj^  =  cosh  (&  -  <£0)/cosh  <£0. 
From  (42)  we  easily  deduce  the  relations, 

(1/t)  sn  (tX  #)  =  sn  (M,  £')/cn  («,  #)  =  sinh  0,, 
en  (tX  £)  =  1/cn  (w,  k'}  =  cosh  <£9 

dn  (iu,  k}  =  dn  (u,  £')/cn  (ti,  #)  =  A0  =  cosh  <£o/cosh  (<k 

......  (46). 

Making  use  of  Jacohi's  well-known  relation  in  the  form 

Z(u,  k1}  =  sn  (u,  *')  dn  (u,  jf)/cn  (M,  k'}  -  vul&KK'}  +  tZ  (iu,  *) 

......  (47) 

we  have  from  (28)  and  (46) 


t-*.)J 


=  sinh  <^0/cosh  (<^  -  &)  -  -n-J?'  (^,  JF)/(2  JOT) 
Noting  that  Legendre's  rektion  gives 

EIK  +  E'/K'-1=$*I(KK")   ...............  (48) 

we  obtain  from  (27) 

E(+,  k'}  =  sinh  ^/cosh  (^  -  ^0)  +  (1  -  E/ST)  F(+t  k'}  +  iZ(iu,  k\ 
or,  finally,  making  use  of  (27),  (42)  and  (44), 
,  k'}  =  sinh  <£0/cosh  (<^  -  «^0) 


—  (<^  sinh  &  +  Cj  sinh  ^>2  +  ...  +cwnili4fe+...)  ...(49). 

(ii)  Calculation  of  the  Jacobian  Functions  sn(u,  k'\  cn(«,  Jf), 
dn  (u,  k'\  etc.  in  terms  of  the  argument  u. 

Compute  the  A.G.M.  scale  (a0=  1,  b0  =  k'}  and  calculate  <frn  from  the 
relation  <f>n  =  2*anu,  n  being  such  that  cn  is  less  than  the  small  quantity 
determining  the  order  of  accuracy  of  the  calculations.  The  hyperbolic 
recurrence  formula 

sinh  (2^._!  -  £,)  =  (eJcO  sinh  <£. 

enables  us  to  calculate  successively  <£„_!  ...  <£8,  &,  and  finally  <£<>>  in 
terms  of  which  we  have  from  (46) 

sn(«,  i')  =  tanh<£o,  en  (u,  £')  =  1/cosh  <^0,   dn  («,  k1}  =  1/cosh  (^  -  «£0) 

......  (50). 

Z(u,  k')  may  then  be  calculated  from  (47)  and  (44). 
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(iii)  Calculation  of  the  Jacobian  Functions  of  imaginary  argument, 
sn  (iu,  k),  en  (iu,  k),  dn  (iu,  k),  etc. 

The  procedure  outlined  above,  together  with  equations  (46),  obviously 
enables  us  to  calculate  the  Jacobian  functions  of  imaginary  argument 
and  modulus  k,  while  formulae  (44)  and  (45)  serve  for  the  functions 
Z(iu,  k)  and  ®(iu,  k). 


SECTION  VI 

GAUSS'  SCALE  OF  INCREASING  AMPLITUDES 
(i)  Circular  Scale. 

In  connection  with  the  A.G.M.  scale,  Gauss  made  use  of  the  recurrence 
formula* 

tan  Xn+i  =  V(««+i/^+i)  tan  2xn   ...............  (51). 

It  is  easily  proved  that  the  above  formula  is  connected  with  Landen's 
scale  (17)  by  the  relation 

tan  Xn  =  *J(bn/an)  tan  fa  .....................  (52) 

or 

fa+l-fa=%Xn-l  ........................  (53). 

In  the  limit,  as  n  tends  to  infinity,  (bnjan}  approaches  unity,  so  that 
Xn  and  fa  both  tend  to  the  same  limit.  It  thus  follows  from  (25)  that 
if  we  form  the  A.G.M.  scale  (a0  =  1,  b0  —  k'\  and  starting  with  x<>  deter- 
mined by  tanxo  =  \/(^o/»o)  tan  <f>,  compute  successively  Xn  /to>  •••  X»  •••  » 
we  have 

^OM)=(i/<Kx»/2«)   ..................  (54). 

From  (52)  and  (53)  it  is  a  simple  matter  to  derive  successively 
fa  ,  fa  .  .  .  fa  .  .  .  for  use  in  the  formulae  of  Section  IV.  When  these  angles 
are  required,  it  is,  however,  simpler  to  make  direct  use  of  Landen's 
recurrence  formula  (17)  or  (18). 

(ii)  Hyperbolic  Scale. 

It  is  evident  that  formulae  (51)  and  (52)  may  be  written  in  the 
hyperbolic  form 

tanh  Xn+i  =  «/(o»+i/6»+i)  tanh  2xn  ...............  (55), 

tanh  xw=.s/(6n/an)  tanh  </>„     ..................  (56). 


*  Gauss,  Get.  Werke  (1866),  Bd.  in.  p.  388.  To  compare  Gauss'  results  with 
those  of  the  present  paper  we  note  that  <£n=2nFn,  xn=2npn+1'  I4  then  apPears 
(p.  392)  that  formula  (26)  of  the  present  book  was  also  obtained  by  Gauss,  inde- 
pendently of  Legendre  (foot-note  on  p.  9). 
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As  n  increases  indefinitely,  x»  and  <£„  tend  to  the  same  limit,  so  that 
if  we  write 

tan  <f>  =  sinh  </>„  or  sin  <f>  =  tanh  <f>0  ...............  (57) 

it  follows  from  (42)  that 

JP(0,*')  =  (l/«-)(^/2")  =  (l/o.)(X^2-)  ............  (58). 

Noting  that,  as  before, 

<k-n-<£»  =  2xn-i    ........................  (59) 

it  is  a  simple  matter  to  derive  successively  <&,  fa,  ...  <£„  ...  for  use  in 
the  formulae  of  Section  V.  In  such  cases,  however,  it  is  simpler  to 
make  direct  use  of  the  hyperbolic  recurrence  formulae  (37)  or  (38). 

(iii)  Calculation  of  1C,  E'  and  the  name  q-e~vK!K. 

If  we  write  <f>  =  ^ir,  it  follows  from  (57)  that  tanh  <t»0—l,  and  from 
(56)  that  tanh  x0  =  VC^oAO-   From  (55)  we  deduce  that 

tanh  2x<>  =  &i/«i>  tanh  2xi  =  k/o,,  ...  tanh2x»  =  fe»+1/aw+1, 
and  thus      2X»  =  \  log  {(a,+1  +  k+O/fat+i  -  6m+1)}  =  \  log  (a,+2/e,,+2) 
or  (x./2»)  =  (£T  8  log  (tWv*)- 

It  follows  from  (58)  that  as  n  is  indefinitely  increased,  (£)"  log  (a^/c,,) 
tends  to  a  finite  limit  given  by 

anF(fr,  F)  =  frK'lK  =  $  log  (I/?)  =  (^^(o./c.)  ...(60). 

Making  use  of  the  relation  (13),  a^Ca^^c2,^,  we  have 

log  (an/Cn)  =  2  log  ^a./a,-,)  +  2  log  (a,.,/^-!), 
which  gives  by  successive  application, 
log  (a,/cn)  = 

2  log  (2an/o,,-i)  +  4  log  (2aw_1/a»_2)  +  .  .  .  +  2"  log  (20,/a,)  +  2*  log  (a,yc0). 

Inverting  the  order  of  the  series,  we  obtain 


+  log  («,/<«,)  +  J  log  (02/Oj)  +  ...  +  (i)"-1  log  (a 
As  n  is  increased  indefinitely  we  obtain  the  series 


(61). 


Making  use  of  Legendre's  relation  (48)  in  the  form 


and  using  (27),  we  obtain  the  following  useful  formulae, 
E'  =  i  (I/a,)  (cc2  +  2c,2  +  4c23  +  .  .  .  +  2"cn2  +  .  .  .  ) 

x  {log  (40,/d)  -  log  (a,/aO  -  £  log  (02/03)  -  -  -  -  }  +  an  ......  (62), 


x  {log  (40,/c,)  -  log  (eh/as)  -  £  log  (aj/Oj)  -  .  .  .  }  -  an  ......  (63) 
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for  the  calculation  of  the  complete  elliptic  integrals  of  complementary 
modulus  in  terms  of  the  A.G.M.  scale  («0  =  1,  b0  =  k'}. 

The  nome  q  is  denned  by  the  relation  log  (Ijq)  =  irK'/K,  from  which 
it  follows  that 

c1)    .........  (64) 


which  is  more  convenient  for  numerical  computation  than  (60).  It  may 
be  noted  that  the  entire  series  of  calculations  may  be  carried  out  with 
ordinary  logarithms  to  the  base  10,  enabling  q  and  powers  of  q  to  be 
evaluated  directly*. 

The  complementary  nome  q'  =  e~vK<K'  is  obviously  related  to  q  by  the 
formula 

lOg  £. 


SECTION  VII 

LANDEN'S  SCALE  OF  DECREASING  AMPLITUDES 

(i)  Calculation  of  the  Elliptic  Integrals  of  the  First  and  Second  kinds, 
modulus  k'. 

If  we  denote 

Are  =  J(an*  cos2  \l/n  +  cj  sin2  $n)  ...............  (65), 

the  recurrence  formula 

tan  (^  -  1/^+1)  =  (cn+1/an+1)  tan  i^n+1  ............  (66) 

may  be  written  in  either  of  the  forms 

sin  (2i/wi  -  <K)  =  (bnfan)  sin  fa,  cos  (2fa+l  -  fa}  -  ATC/ara  .  .  .  (67). 
From  these  we  deduce 

Aw  +  bn  cos  fa  =  2AW+1,  AB  -  bn  cos  fa  =  2«n+1cn+1/An+1  ...(68), 
which  give  by  differentiation 

%/A0  =  ^1/A1  =  ^2/A2=  •••  =  dfa/*n-=  .........  (69). 

If  we  construct  the  A.G.M.  scale  (a0  =  l,  b0-k')  we  notice  that  as  n 
increases,  tyn  tends  rapidly  to  a  finite  limit.  If  we  start  with  i/^0  =  $, 
A0  =  A  (\j/,  k'},  we  have,  since  cn  ultimately  tends  to  zero, 

*  The  limit  derived  in  (60)  does  not  appear  to  be  generally  known,  although  an 
equivalent  formula  is  given  by  Legendre  (Fonctions  Elliptiques,  t.  i.  pp.  94,  101). 
Although  formula  (60)  does  not  appear  to  have  been  explicitly  stated  by  Gauss, 
the  series  (61)  is  found  among  a  number  of  formulae  given  in  his  posthumous 
papers  (Gauss,  Ges.  Werke,  Bd.  in.  p.  377).  Neither  (60)  nor  (61)  appears  to  have 
been  noticed  by  Jacob!  or  by  subsequent  writers  on  elliptic  functions. 
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,ir)  =  (l/an)  / 
Jo 

a  useful  formula  due  to  Legendre*. 
From  (68)  we  easily  obtain  the  relation 

A»s  +  ft.  A,  cos  *n  =  2A 
which  gives 


Writing  down  the  system  of  equations  for  n  =  0,  1,  2,  ...  n,  multiplying 
successively  by  1,  2,  4,  ...  2",  adding  and  integrating,  we  obtain  the 

result 


ro  /- 

/     A0<fy0-2»+1 
/o  .'o  Jo 

-(&„  8^0  +  2^8111^1+...  +2"6,sin^,)  ......  (71). 

Adding  to  each  side  the  identity 


/•K+l 

2"+1  bn+l  cos  ^d^  =  6»+1  sin  ^w+1  +  (2"+I  -  1)  bn^  sin  *»+, 

Jo 

and  re-arranging  the  terms,  (71)  may  be  written 

rto  rt*+i 

\     A0^0-2'+1/        (A,+1-6n+1  cos  ^O^+i 
Jo  Jo 

rto 
=  %  (c0a  +  2Ci2  +  .  .  .  +  2wcn2)  I     cfy0/A0  +  6W+1  sin  <jfn+1  +  (^  sin  ^  -  b0  sin  $») 

,  sin  ^2  -  6X  sin  ^)  +  .  .  .  +  (2"+1  -  1)  (bn+l  sin  <ffn+1  -  bn  sin  ^) 

......  (72). 


It  isjwf  difficult  to  establish  the  identity 
Vi  sin  tn+i  -  bn  sin  «^n  =  2e,,+2  tan  (2^,+a  -  </w)/cos  (2^m+3  -  ^n+2) 

......  (73). 

As  n  increases  indefinitely,  an+1  and  bn+1  tend  to  the  same  limit  an, 
while  cn+1  tends  to  zero.   It  follows  that  the  second  integral  on  the  left- 
hand  side  of  (72)  tends  to  zero,  so  that  we  finally  obtain,  making  use 
of  (70), 
Efy,  k")  =  an  sin  ^  +  |  (c02  +  2Cl2  +  ...  +  2"c,a  +  ...) 

x  (I/a,)  log  tan  (£*•  +  |>H)  +  2  (^.,  c.)  ......  (74), 

where  for  the  sake  of  brevity  we  denote 


(75), 


2^    c^-Zc  -        ,  6c  -        . 

2C2cos(2^3-vA2)        3cos2,A4-^ 


*  Legendre,  Functions  Elliptiques,  t.  r.  chap.  xrx. 
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a  useful  formula  for  the  numerical  calculation  of  the  second  elliptic 
integral*. 

Making  use  of  (27),  (28),  (48)  and  (70),  we  obtain  from  (74)  the 
following  interesting  expressions 

Zty,  k1}  +  I*  Fty,  k'}l(KK'}  =  an  sin  ^  +  S  (^,  cn)  ...(76), 

2 

or,  in  terms  of  u,  by  (70)  we  have 

sin  \f/n  =  tanh  (anu),  \ 

and  Z  (u,  k'}  +  anujK'  =  an  tanh  (anu)  +  2  («//„,  c»)  |      •  •  •  (7  ?)• 

2  J 

Jacobi's  relation  (47),  taken  with  equation  (76),  gives 

(1/0  Z  (iu,  k}  =  tan  $0  cos  (2^  -  «/rc)  -  an  sin  ^n  -  2  (</^,  cn)  (78). 

2 

The  addition  theorem 


Z(u)  +  Z(v)-Z(u  +  v)  =  tf  sn  (M  +  v)  sn  M  sn  v  ......  (79), 

with  iu  for  u  and  t>  =  ±  (K+  iK'}  gives  two  interesting  formulae 
iZ(iu  +  K+  iK',  K)  =  an(l+  sin  \l/n) 

...(80). 


If  we  notice  that 

I     sin  \f/nd\j/nl&n  =  (l/bn)  log  (an+\l  &n+\)  =  i 
Jo 

the  result  of  integrating  (72)  written  in  the  form 

Z(u,  k'}  +  anu/K'  =  an  sin  \}/n  +  2  (2ra+1—  1)  (6M+1  sin  ^n+1  -  6W  sin 

o 

gives  rise  to  a  formula  for  the  numerical  calculation  of  ®  (u,  k'\ 


in  which,  if  desirable,  we  may  write  sec  \f/n  =  cosh  (anw). 

(ii)  Calculation  of  the  Jacobian  functions  sn  (M,  ^'),  CU(M,  £'),  dn  (w,  ^'), 
e#c.  in  terms  of  the  argument  u. 

Form  the  A.G.M.  scale  (a0=l,  bQ  =  k'}  and  determine  i^re  from 
sin  ^n  =  tanh  (anu),  n  being  such  that  cn  is  less  than  the  small  quantity 
determining  the  order  of  accuracy  of  the  calculations. 

*  The  series  (71)  was  obtained  by  Legendre  (Fonctiom  Elliptiques,  1. 1.  chap,  xxi.) 
and  left  by  him  in  this  non-convergent  form.  As  far  as  the  writer  is  aware,  the 
convergent  series  (75)  is  new. 
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From  the  recurrence  formula  (66),  tan  (i/rn_!  -  i/rn)  =  (cn/an)  tan  $*, 
compute  i/^-i ...  ifrz,  ^1}  and  finally  i^0,  in  terms  of  which 
sn  («,£')  =  sin  i/^,     cn(«,  k')  =  cos$0,     dn (u,  k'}  =  A0  =  cos (2^ - ^0) 

(82); 

')  may  then  be  calculated  from  (77)  and  ®(u,  k'}  from  (81). 


SECTION  VIII 

THE  HYPEEBOLIC  SCALE  OF  DECREASING  AMPLITUDES 

(i)  Calculation  of  the  Elliptic  Integrals  of  the  First  and  Second  kinds, 
modulus  k. 

If  we  write  i*fr  for  ty  throughout  the  formulae  (65)  to  (69),  (65)  becomes 
A.=  V(a.fcoshV.-c.1sinha^ll)  ...............  (83), 

while  the  recurrence  formula  (66)  may  be  written 

tanh  (</rn  -  ^.+1)  =  ((VnMi+i)  tanh  ^,+1     .........  (84), 

or,  more  conveniently, 

sinh  (2^n+1  -  ^B)  -  (&B/O  sinh  ^n  ,     cosh  (2tin+l  -  ^B)  =  Aw/a»  .  .  .  (85). 
As  in  Section  VII  (i),  it  is  readily  proved  that 

<fyo/A0  =  ^I/A!  =  rffc/A,  -  ...  -  ^n/An  =  .........  (86). 

If  we  construct  the  A.G.M.  scale  (a0  =  1,  b0  =  k'}  and  write 

tan  $  -  sinh  ^0     ........................  (87), 

it  is  readily  proved  that 

A0  =  A  (</r,  &)/cos  ^  and  <fy0/A0  =  <^/A  (^,  k}    ......  (88). 

Thus,  writing  u  =  F(ty,  k\  we  have 

sn  (u,  k)  =  sin  ^  =  tanh  ^0,     en  (u,  k}  =  cost  =  sech  i/O 
dn  («,  *)  =  A  (^,  *)  =  A0/cosh  ^0  =  cosh  (2^  -  ^0)/cosh  ^J"    "'  ^89^ 
As  n  increases,  the  ^'s  tend  to  a  limit  \f/n,  while  cn  tends  to  zero. 
Thus  the  result  of  integrating  (86)  gives 


1  (sinh  ^»)    ...(90). 


If  we  notice  that 
r°A0^0=f  A( 

70  JO 

-  r  dn2  («,  *)  dujctf  (u,  &)  -  ["  dn2  (i  a,  >t')  rfw, 
/o  ./o 
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we  have  from  (29) 

/    A0  c?ft,  =  ( 1  /O  Z  (iu,  &')  +  (E'lK'}  u. 
Jo 

It  is  easily  seen  that  equation  (72)  holds  good  when  ityn  is  written 
for  ft,.   Hence  if  we  write 

•$'  /          \  —  f>    tanh  (2fts  ~  ft)  ,  e  „  tanh  (2^3  -  r'- N 
2  WM^)"   Ca  Cosh  (2^3-^ 


a     ,  , 

cosn  ^yn+3  — 

and  make  use  of  (47)  and  (48),  we  obtain  after  some  reductions 
Z  (f,  k)  =  tanh  ft,  cosh  (2ft  -  ft,)  -  an  sinh  ft,  -  2'  (</rn,  cn)  . .  .(92). 

2 

The  relations  (27)  and  (28)  enable  us  to  calculate  the  second  elliptic 
integral  Ety,  k\ 

(ii)  Calculation  of  the  Jacobian  Functions  sn  (u,  k),  en  (u,  fc),  dn  (u,  &), 
etc.  in  terms  of  the  argument  u. 

Compute  the  A.G.M.  scale  a0  =l,b0  =  k'  and  determine  f»  from  (90), 
which  may  be  written 

sinh  \l/n  =  tan  (anu)     (93), 

n  being  such  that  cn  is  less  than  the  small  quantity  determining  the 
order  of  accuracy  of  the  calculations. 

We  then  calculate  the  amplitudes  i/'n+1  ...  i^2,  «Ai,  "Ao  from  the  recur- 
rence formula 

tanh  (i/^-!  -  ij/n)  =  (cn/Ofc)  tanh  ^n. 

Equation  (89)  then  gives 

sn  (u,  k)  -  tanh  \j/0 ,    cn  (u,  K)  =  sech  \l/0] 


(94). 

dn  (u,  k)  =  cosh  (2^  -  ^0)/cosh  \l/0 


Equations  (91)  and  (92)  then  enable  us  to  calculate  ^(w,  A;).    To 
compute  ®  (w,  A;),  we  note  that  (92)  may  also  be  written 
Z(u,  k}  -  sn  (u,  k~)  dn  (u,  k)/cn  (u,  k)  —  an  sinh  ftt 

-  2  (2"+1  -  1)  (&n+1  sinh  ^w+1  -  bn  sinh  ftB). 

o 

Integrating  this  expression  with  respect  to  u  and  noticing  that 

sinh  ftn^n/Aw  =  (!/&„)  log  cosh  (2ft,+2  -  ^+1), 
/o 

we  obtain  the  following  series  for  ®  (u,  k), 

^cosl 


if 

/i 


Inr  '         -  loir 

10g  g 


«+1  - 


- 
cosh  f  cos 

......  (95). 
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SECTION  IX 

ON  THE  NUMERICAL  COMPUTATION  OF  THE 
THIRD  ELLIPTIC  INTEGRAL 

As  defined  by  Legendre,  the  third  elliptic  integral  may  be  written 


ns(n,  £,<£)=         -  .         A  /  .    I.N  ............  (96). 

Jo  (l+»sin2</>)A(<£,  k) 

The  parameter  n  is  taken  as  real,  and  may  take  any  value,  positive 
or  negative. 

Writing  u  =  jF(<£,  k),  n=-k*  sna  (a,  k\  Jacobi  considers  the  integral 

,„„. 
...(97). 

v     ' 


o  -s 

Legendre's  integral  is  connected  with  Jacobi's  by  the  relation 

H3  (»,  k,  <f>)  =  M  +  sn  a/(cn  a  dn  a)  .  n  (w,  a)    ......  (98), 

and  in  this  book  is  differentiated  from  it  by  the  suffix  3. 

In  order  to  include  all  real  values  of  n  in  (96),  a  is  not  restricted  to 
be  real,  and  the  treatment  is  usually  divided  into  four  cases  : 


Case  I.  n  negative  between  0  and  - 
Case  II.  n  negative  between  -  1  and 
Case  III.  n  negative  between  -tr  and  —  1 


ase       .   n  negatve    etween  -tr  and  —  11    ~.      7     r 
Case  IV.   n  positive  between  0  and  +  oo   J  ( 

CASE  I.   n  negative,  between  0  and  -A2.   Hyperbolic  case. 
We  write    n  =  -  Tc1  sns  (a,  k}  =  -  tf  sin2  6,  where  0  <  a  <  K  ...(99). 
From  (97)  and  (98)  we  have 


1  __  SD(^^)  _  w  r®( 

*  Cn  (a,  *)  dn  (a,  *)  10g  L®  ( 


__  _ 

Cn  (a,  *)  dn  (a,  *)  ®  (a  +  «, 

Compute  the  A.G.M.  scale  (a0  =  1,  b0=  k'),  and  starting  with  <f>0  =  0,  make 
use  of  the  recurrence  formula  (17),  tan  (^>n+1  -<^w)  =  (bjaj  tan  <£„.   We 
then  have  sn  (a,£)=sin  <f>0,  en  (a,  k}  =  cos  <£0>  dn(a,  k)  =  cos  <£o/cos(<k-<£0), 
and  Z^a,  ^)  =  d  sin  <&  +  c2  sin  ^>2  +  .  .  .  +  cn  sin  <£„  +  .  .  .  ,  while  ana  = 
We  thus  have 


C 


1  + 

log 


®  ( 

2—2 
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The  calculation  of  u  by  the  recurrence  formula  (17),  commencing  with 
<Ao  =  <£>  presents  no  difficulty.  The  most  convenient  method  of  computing 
the  last  term  of  (101)  is  to  make  use  of  Jacobi's  series 

©(#)=  l-2qcos2anx  +  2q*cos4:anx-2q9cos6ana;  +  ......  (102) 

in  which  we  have  written  an  =  \^IK.   The  nome  q  is  most  conveniently 

00 

calculated  from  (64),  log  q  =  2  (^T~l  l°g  («n/«n+i)  -  log  (^/c,). 

If  we  write  <£  =  \tr,  the  last  term  of  (101)  vanishes,  and  we  obtain  for 
the  complete  elliptic  integral  of  the  third  kind  the  expression 

II3  (n,  k,  |TT)  =  K\l  +  tan  <£„  sec  <f>0  cos  (<k  -  <£<>)  fci  sin  <h  +  c2  sin  <£2  + 

...+cn  sin  <£„+...}]  ......  (103). 

We  may  also  make  use  of  the  hyperbolic  scale  of  decreasing  amplitudes 
and  the  formulae  of  §  8.  The  final  results  are  not,  however,  so  con- 
venient for  computation  as  (101)  and  (103). 

CASE  II.   n  negative,  between  —  1  and  —  <x>  .   Hyperbolic  case. 

We  write 

n  =  -  A2  sn2  (a  +  iK',  A)  =  -  1/sn2  (a,  k)  =  -  I/sin2  6,  where  0  <  a  <  K 

......  (104). 

Equation  (98)  becomes 


=u  --        _n(«,a  +  fZ?)  ...(105). 


2         ,         --        _ 
o  sn2a-sn2w          en  a  dn  a 

We  have 

Z(a  +  iK')  =  Z(a')  +  cnadna/sna-^7ri/K  ......  (106) 

and  ®(z  +  iK')  =  iBH(z)  .....................  (107), 

where  \^B  =  \^K'IK-\^izlK. 

H(z)  is  Jacobi's  Theta-function  given  by  the  expansion* 

H(z)  =  2gi  (sin  anz  -  (f  sin  3anz  +  q6  sin  5anz  -  q™  sin  lanz  +  ...) 

......  (108). 

Writing  successively  (a  -  u}  and  (a  +  u]  for  z  in  (107),  and  a  +  iK  for  a 
in  (97),  we  obtain 

H  (u,  a.  +  iK')  =  uZ  (a,  k}  +  U  Cn  a  dn  a/sn  a  +  |  log  [H  (a  -  tl)IH(a  +  tt)], 

and  hence  from  (105),  u  <  a, 

H3  («,  k,<f>)  =  -  sn  o/(cn  a  dn  a)  [uZ  (a,  A)  +  £  log  {H  (a  -  u)/H  (a  +  «)}] 

......  (109). 

*  Jacobi,  Fundamenta  Nova,  §§  62,  63  ;  Ges.  Werke,  1881,  Bd.  n.  pp.  229-231  ; 
Cayley,  Elliptic  Functions,  2nd  Edition  (1895),  Chap.  vi. 
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For  numerical  calculation,  form  the  A.G.M.  scale  (a0  =  l,  b^  =  1i},  and 
starting  with  <f>0  =  0,  make  use  of  the  recurrence  formula  (17), 

tan  (<k+i  -  <k)  =  (&»K)  tan  <f>n  . 
We  then  have  ana  =  <f>*j2*,  and  as  in  Case  I,  we  derive  finally 

n,(n,  *,<*>) 

-  -  u  tan  <£0  sec  <f>0  cos  (<&  -<£<>)  (ci  sin  </>!  +  ^  sin  <k  +•••+  Cm  sin  <£„+...} 
-  ?  tan  <£„  sec  <f>0  cos  (&  -  <£„)  log  [H(a  -  ii)/H(a  +  uj]   ......  (HO). 

Starting  with  <£„  =  <f>  and  making  use  of  the  same  recurrence  formula  as 
above,  the  computation  of  u  presents  no  difficulty.  As  before,  the 
calculation  of  the  last  term  of  (110)  is  best  achieved  by  computing  the 
nome  q  by  (64)  and  making  use  of  (108). 

If  we  write  <£  =  |ir,  the  last  term  of  (110)  vanishes,  and  we  obtain  for 
the  complete  elliptic  integral  of  the  third  kind  the  expression 


=  -  JTtan  <j>o  sec  <£0  cos  (<&-  <£0)  fa  sin  &  +  c^  sin  <k  +  ... 

......(111). 

When  M>a,  it  will  be  noticed  that  the  integral  (105)  has  a  singularity  at 
tt=o,  where  theintegrandbecomes  infinite.  Atthe  same  time  the  logarithm 
of  a  negative  quantity  appears  in  (110),  since  H(a  —  'u)  =  —  H(u  —  a). 
It  is  necessary  in  this  instance  to  define  the  path  of  integration,  re- 
garding (105)  as  an  algebraic  integral  obtained  by  writing  z  =  sn  u.  Near 
the  point  z  —  sn  a  this  integral  is  approximately 

n3  =  sn2  a/(cn  a  dn  a)  .  /  <&/(sn2  a  -  2s), 

which  integrated  round  a  small  semi-circle  below  the  real  axis  gives  rise 
to  the  imaginary  term  —  \tri  sn  a/(cn  a  dn  a).  If  we  define  the  path  of 
integration  of  the  algebraic  form  of  the  integral  as  consisting  of  this 
small  semi-circle  and  the  remaining  portions  of  the  real  axis,  the  above 
imaginary  term  cancels  the  term  —  sn  a/(cn  a  dn  a)  .  ^  log  (—  1)  appearing 
on  the  right-hand  side  of  (110)  when  we  write  H(a  —  tt)  =  —  H(u-a). 
"With  this  meaning  attached  to  the  path  of  integration,  we  may  write 
the  incomplete  integral  in  the  form,  u  >  a, 
Hs  (»,  £,<£)  =  -  Sn  a/(cn  a  .  dn  a)  [uZ  (a,  *)  +  i  log  {H  (u  -  a)/H(u  +  a)}] 

......  (112). 

Numerical  computation  proceeds  in  the  manner  already  described. 

CASE  III.   n  negative,  between  -A?  and  -  1.    Circular  case. 

First  Method. 
We  write 

«  =  -Fsn2(ia  +  Jg;^)  =  -A2/(ln2(a,/fc')  =  -/ta/(l-/t'3sin2^)    ...(113), 
whence  sin  6  =  sn  (a,  k'\  0  <  a  <  K1  . 
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Writing  ia  +  K  for  a  in  (98),  we  find 

H3  (n,  k,tf  =  u  +  i  dn  (a,  #)/[#*  sn  (a,  £')  en  (a,  #)]  .  n  (u,  ia  +  IT) 

......  (114), 

while  (97)  gives 


H  (u,  ia+K}  =  uZ(ia  +  -£")  +  £  log  [®  (t'a  -  w  +  Jf)/0  ^a,  +  w  +  JT)] 

......  (115). 

Making  use  of  the  addition  formula  (79)  for  Z(u),  we  find  after  some 
reductions 

iZ  (ia.  +  XT)  =  iZ  (ia,  k}  +  V  sn  (a,  #)/[cn  (a,  V)  dn  (a,  #)]  ...  (1  1  6), 
or,  from  (47),  in  terms  of  a  real  argument 

iZ(ia  +  K}  =  Z(a,  k'}  +  1  wo/(JCr')  -  A;'2  sn  (a,  #)  en  (a,  #)/dn  (a,  #) 

......  (117). 

Introducing  Jacobi's  Theta-function  ®!  (z)  defined  by  the  relations 

®(z  +  K}  =  ®l(z)    .....................  (118), 

®!  (z)  =  l+  2q  cos  2anz  +  2cf  cos  4are2;  +  2/cos  6anz  +  ......  (119), 

we  may  write  (114)  in  either  of  the  forms 

-          dn  (a'  ^ 

-      -  - 


(*,  I  ^  dn  (a,  k') 

ls  (  »,  K,  <p)  —  775  -  -  -     - 

2 


75        -,  -  Yi\  -  7       77\ 

k  2  sn  (a,  k)  .  en  (a,  k) 


For  numerical  computation  we  may  make  use  of  two  sets  of  formulae, 
according  as  we  employ  circular  or  hyperbolic  functions. 

(i)  Use  of  circular  recurrence  formulae. 

Form  the  A.G.M.  scale  (a0  =  1,  b0  =  k').  Commencing  with  </>0  =  <£,  make 
use  of  the  recurrence  formula  (17),  tan  (<£„+!  -  <£n)  -  (bn/a>n)  tan  <£„,  and 
thus  derive  anu  =  <£n/2n  =  $. 

Commencing  with  ^0  =  ^  employ  the  recurrence  formula  (67), 
sin(2«/'B+1-^)  =  (6w/an)sm^n,  giving  «wa  =  logtan(^7r  +  |-fn).  The 
application  of  (47)  and  (76)  to  both  (120)  and  (121)  gives  us  the  formula 

TT  /     z.  J.N    2  cos  (2^  -  i/^o)  r    .       .  -/,M     1-1     ®i(«-t 


......  (122) 

for  the  numerical  calculation  of  the  incomplete  integral,  and  the  formula 
H8  (w,  k,  Jir)  =  2  (  JT/£'2)  cos  (2^ 
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The  peculiar  difficulty  arising  in  the  numerical  computation  of  the 
circular  cases  of  the  incomplete  elliptic  integral  of  the  third  kind  lies  in 
the  evaluation  of  the  Theta-functions  of  complex  argument.  Writing 
z  =  u  +ia  in  Jacobi's  expansion  (119),  the  series  may  be  resolved  into 
real  and  imaginary  parts*.  If  we  write 

2g  sin  2«n«  .  sinh  2a<ta  +  2<y*  sin  4an«  .  sinh  4a»a  +  .  .  . 
~  1  +  2q  cos  2anw  .  cosh  2a»a  +  2q*  cos  4a»tt  .  cosh  4ana  +  .  .  . 

......  (124), 

we  easily  find  log  [®i  (u  -  to)/®i  (»  +  •°)]  =  2**>  fr°m  which  it  follows 
that  (122)  may  be  written 

n3  (n,  k,  </>)  =  (2/*'a)  cos  (2«Ai  -  ^0)/sin  2^0  [«  K  sin  ^,  +  1  ty)}  -  *] 

2-  .....  (125). 

Powers  of  the  nome  q  arising  in  (124)  are  most  conveniently  calculated 
from  (64), 

00 

log  q  =  2  (I)"'1  log  (o»/a»+1)  - 


The  use  of  hyperbolic  functions  in  (124)  may  be  avoided  by  noting 
that  if  we  write 

<£  =  02*."  =  tan2  (i»  +  itfr.)  =  (1  +  sin  ^)/(l  ~  s^  fc,), 

Q2+Q-2)        2cosh4ana=  Q*+  Q-^ 
Q2-Q-2j  '      2sinh4ana=  (f-Q-*)  '  " 

from  which  it  follows  that  (124)  may  be  written 
^ 


......  (124  Us}. 

(ii)  Use  of  hyperbolic  recurrence  formulae. 

Write         »  =  -  #  sn2  (ta  +  ^,  *)  =  -  A»/(l  -  A;'2  tanh2  ^O  ......  (126), 


whence  sn  (a,  k'}  =  tanh  ^>D',  0  <  a  <  K'. 

From  the  hyperbolic  recurrence  formula  (37) 

tanh  (^'n+1  -  <£„')  =  (bn/an)  tanh  ^ 
we  have 

an  a  =  <A»72"  =  *',     sn  (a,  £')  =  tanh 

en  (a,  K}  =  sech  ^>0'  ,    dn  (a,  k'}  =  sec 


*  This  procedure  is  dae  to  Legendre,  Fanetiant  Elliptiques,  Deaxieme  Supplement 
(1828),  §u.  p.  146  ef  teq. 


24  NUMERICAL  COMPUTATION  OF  THIRD  INTEGRAL 

Thus  from  (44)  and  (120)  we  obtain 

TT   (      Z-   £\—   1   _  COsh  <j>0'  _ 

13  (n'  *  V  - 


_  _ 

F2tanh</.0'cosh(<k'-4»0') 


f    wsinh^o'  ,£       •  i    ,  /     1-1     ®i(u~ia)~\ 

X      •—-.,,,        ,  ,x  -1*2  cn  Sinh  <t>n  +£«  log  =*-)  -  T-C         ...(127). 

Lcosh  (&'-<£„')      T  6®1(M  +  za)J 


To  calculate  u  we  make  use  of  the  recurrence  formula  (17) 

tan  (<£„+!  -  <£„)  =  (6«/o»)  tan  <£„,  giving  «nw  -  4>n/2w  =  $. 

Equation  (124)  then  becomes 

,          _    2q  sin  2$  sinh  2$'  +  2^  sin  4$  sinh  4<I>'  +  2qg  sin  6$  sinh  6$'  +  .  . 
1  +  2q  cos  2$  cosh  2$'  +  2<?4cos  4$  cosh  43>'  +  2^9  cos  64>  cosh  6$'+  . 

......  (128). 

Equation  (127)  then  gives 

rr  f     z-  rk\  -  J:  cosh  ^>0' 

3  ^  *»  ^;  " 


and 

2         k'~  tanh  <£„'  cosh  (<£/  —  </>0')  Lcosh  (<£/  —  $<>')     i 

(129  to). 

Second  Method. 
If  we  write  ft  =  K'  -  a,  (113)  takes  the  form 


(130). 

We  easily  find 

n-  /.  i  OA  _  dn  (A  #) 


j  ,•  log 

B 


sn  cn 


.„,     ._N     TT       .  .,         l  /ioiN 

1^  -  1^)  -  —  +  1  1  log  ^  ^  _  .          .........  (131) 

dn  08,10 


.       . 

-(  }i 


where  ^Ti  (2;)  is  Jacobi's  Theta-function  defined  by  the  series 

HI  (z)  =  2#i  (cos  anz  +  q2  cos  3anz  +  <f  cos  5anz  +  g-12  cos  7an«  +  .  .  .) 

......  (133). 
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Writing  *'  =  |«log[J71(«  +  t/8)/J5r1(«-ij8)]      .........  (134), 

we  find 

,       .,  _  sin  anMsinha,t/8+y8sin3anMsinh  3anP  +  qBsm5anusinh5anft+  ... 
cos  an  M  cosh  anf$+q*  cos  3an  M  cosh  3anft+<f  cos  5a»  u  cosh  5an  fi+... 

......  (135). 

The  numerical  computation  may  be  carried  out  by  two  methods. 

(i)   Use  of  circular  recurrence  form  ulae. 

Starting  with  \j/0  =  &  in  (130),  make  use  of  the  recurrence  formula  (67), 
sin  (2<A»+i  -  &•)  =  (&»/«»)  sin  &»•  From  (80)  we  have 

iZ  (K+  iK'  -  tf)  =  a.  (1  -  sin  0.)  -  2  (f  ), 
while 

sn  (ft  k'}  =  sin  ft,,  en  (&  A')  =  cos  i^0  ,  dn  (^  *')  =  cos  (2fc  -  ^0). 
Thus  (131)  becomes 

n.(«^^)  =  «  +  2^^-^[^-t,{o.8infc+aW}]    ...(136). 

As  before,  «  is  most  conveniently  calculated  from  the  recurrence 
formula  tan  (<£n+1  -  <£„)  =  (6,/a,)  tan  ^n,  giving  a.M  =  <k,/2"  =  $. 

The  use  of  hyperbolic  functions  in  (135)  may  be  avoided  by  writing 
log  q  =  anfi  =  log  tan  (£*•  +  Jft,),  so  that 

2  cosh  anp=Q+  tf-i.         2  cosh  3«M/?  =  ^  +  tf-«i 

2  sinh  anp  =  q-  q-*}  '      2  sinh  3aw/8  =  Q3  -  Q-*}  '  6 

Thus  (135)  may  be  written 


......  (138). 

Writing  u  =  JTin  (136),  **  =  fa  and  we  have,  since  $v  =  anK, 


(ii)    Use  of  hyperbolic  recurrence  formulae. 

Write  sin  ^  =  sn  Q8,  ^')  =  tanh  <£„',  and  (130)  becomes 

»  =  -  #  sn2  (IT  +  iK'  -  ip)  -  -  1  +  **  tanh2  ft    .  .  .(140). 
Make  use  of  the  hyperbolic  recurrence  formula  (37) 

tanh  (</>'w+1  -  £/)  =  (6,/a.)  tan 
and  we  have 


en  (/3,  A')  =  sech  «£„',    dn  (ft  ^)  =  sec 
Also    -  z'Z  (ift,  K)  =  Cj  sinh  ^'  +  c2  sinh  ^a'  +  .  .  .  +  c,  sinh  £,'+..., 
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so  that  (132)  gives 
H3  (n,  k,  <£) 

1  cosh  <£„' 


k'2  cosh  (&  -$<>')  tanh  </>0' 

x  |~-  uk2  sinh  <f>0'  cosh  (&'-  <£<>')  +  «  2  cre  sinh<£n'  +  *'  |  ...(141). 

As  before,  we  calculate  u  from  the  recurrence  formula 

tan  (<£n+1  -  <f>n)  =  (bn/an)  tan  <£n, 
giving  anu  =  ^n/271  =  <J>.    Equation  (135)  gives  for  *' 

,  _  sinh  <&'  sin  $  +  q*  sinh  3$'  sin  3<E  +  g6  sinh  5$'  sin  5$  -t-  .  .  . 
cosh  3>'cos  $  +  <fcosh  3$'  cos  3$  +  ^8cosh  5$'  cos  5$  +  .  .  . 

......  (142). 

Writing  u  =  K  in  (141),  we  have  *'  =  JTT,  and  since  |TT  =  anK, 

H3  (n,  k,  |ir) 


......  (143). 

CASE  IV.    n  positive,  between  0  awrf  <x>  .    Circular  case. 

First  Method. 

Write    n  =  -  k2  sn2  (to,  A)  =  *»  sn2  (a,  ^')/cn2  (a,  k'}  =  1?  tan2  0.  .  .  (144), 
giving  sin  0  =  sn  (a,  k'},     0  <  a  <  K'. 

Thus 

TT  /      7    j\  Sn  (a,  k'}  en  (a,  &')  P  .     „  .  .  N      ,.,       ® 

n3  („,  k,  *)  =  «  +  -  ^^-  -;  [*^W  +  i«  log 

......  (145). 

®  («)  is  Jacobi's  Theta-function  defined  by  the  series 

©  (z)  =  1  -  2q  cos  2anz  +  2q*  cos  ±anz  -  2#9  cos  Banz  +  ...  (146). 
Writing  V  =  %ilog[®(u-ia)/®(u  +  ia)]     ............  (147), 

we  find 

.  ____  2q  sin  2anu  sinh  2ana  -  2q*  sin  4«MM  sinh  4ana  +  .  .  . 

1  -  2q  cos  2an  u  cosh  2an  a  +  2qf*  cos  4a»  M  cosh  4an  a  -  ... 

......  (148). 

The  numerical  computation  may  be  carried  out  by  two  methods. 

(i)    Use  of  circular  recurrence  formulae. 

Make  use  of  the  recurrence  formula  sin(2i/fn+1-^n)  =  (6n/an)sini/' 
commencing  with  ^  =  0  in  (144).    We  then  have 

sn  (a,  k'}  -  sin  \l/0,     en  (a,  k')  =  cos  ^0,     dn  (a,  k'}  =•  cos  (2^i  -  1/^), 
and  by  (78) 

iZ(ia,  Jc}  =  an  sin  jn  +  2  («A)  -  sn  (a,  #)  dn  (a,  ^)/cn  (a,  A')  ...(149). 
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Hence  (145)  may  be  written 


i 

......  (150). 

As  before,  t*  is  calculated  from  the  recurrence  formula 

tan  (<£n+1  -  <£,,)  =  (&«/«„)  tan  £,, 
giving  o»  M  =  <£n/2n  =  4>,  while  a»a  =  log  tan  (£  ^  +  £  ^n).    We  have 


awa=  .       . 

'    2sinh4a.a=  ^-Q-«J  ' 
where  Q  =  tan  (J  IT  +  i/rw), 

so  that  (148)  becomes 


......  (152). 

Writing  w  =  ^in  (150),  *  =  0,  and 


.--^         t  +  °"  sm  *• 

......  (153). 

(ii)    Use  of  hyperbolic  recurrence  formulae. 

Write  sin0  =  sn(a,£')  =  tanh<£0'  ..................  (154), 

so  that  w  =  -A3sn2(«o,^)  =  >t2tan2e  =  Arsinh2^0'. 

Starting  with  <&,'  given  by  (154),  make  use  of  the  recurrence  formula 

tanhC^m  -^,')  =  (6wK)tanh^'  ............  (155). 

We  then  have 

o»a  =  $.72"  =  ^»',         sn  (a,  >fc')  =  tanh  <&>', 
en  (a,  1Q  =  sech  ft,     dn  (a,  A:')  =  sec 

00 

and  -  t'Z  (to,  K)  =  X  CB  sinh  <£„'. 

Equation  (145)  may  then  be  written 


As  before,  M  is  calculated  from  the  recurrence  formula 

tan  ($m+1  -  $„)  =  (fen/a,,)  tan  «^,, 
giving  Onit  =  $»/2'*  =  <t>.   Equation  (148)  then  gives  for  * 

2<y  sinh  23>;  sin  2$  -  2<y*  sinh  4^'  sin  4*  +  2<y9  sinh  63>'  sin  6O  -  .  .  . 
~ 


......  (157). 
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Writing  u  =  K  in  (156),  *  =  0,  and 

<>'  cosh  (<£/  -  <£o')  "        .  ,   ,  ,~l     /,  K0\ 
—  '2cnsmh</>n'    ...(158). 
cosn  <PO  i  j 


Second  Method. 
If  we  write  fi  =  K'  —  a.  in  (145),  we  have 

n  =  -  tf  sn2  (iK  -  ip)  =  cna  (0,  £')/sn2  (ft  #)  =  cot2  d  .  .  .(159), 
giving  sn  (ft  £')  =  sin  6,     0  <  0  <  JT. 

We  easily  find 


,4Qen(fl,JO  t.  _,.._,    . 

;0  7\  iuZ(iK  -i 

n  (/?,  ^  )         L 


where  ^(«)  is  Jacobi's  Theta-f  unction  defined  by 
H(z)  =  2q*  (sin  anz-qz  sin  Bo^z  +  5^  sin  5anz  -  q12  si 

,,  .4 

It  we  write 


]     .........  (162), 

we  find 

,  _  ,  _  cosareMsinhgn/?  -  <f  cos  3anu  sinh  3an(3  +  g^cos  5anu  sinh  5an/8-  .  .  . 
sin  «nw  cosh  an(3  -  <f  sin  3anw  cosh  3an(3  +  <f  sin  5anu  cosh  5an)8  -  ... 

......  (163). 

The  numerical  computation  can  be  carried  out  by  two  methods. 

(i)    Use  of  circular  recurrence  formulae. 

Commencing  with  ^0=  0in(159),make  useof  therecurrenceformula(67) 

sin  (2^n+1  -  i/O  =  (&„/«„)  sin  ^. 
From  (78)  we  have 

-  iZ(ift,  *)  =  tan  ^0  cos  (2^  -  ft,)  -  «„  sin  ^re  -  2  (ft)  .  .  .(164), 

2 

while  sn(/2,  A')  =  sinft0,    en  (/8,  £')  =  cos  ft0  ,    dn  ((3,  k'}  =  cos  (2ft!  -  ft,)- 
Equation  (160)  then  gives 


......  (165). 

As  before,  u  is  calculated  from  the  recurrence  formula 

tan  (<f>n+1  -  <£„)  =  (6n/an)  tan  <#>„, 
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giving  anu  =  <£n/2"  =  «fr,  while  a*p  =  log  tan  (£TT  +  |i/y).    Thus 
2  cosh  an/3  =  Q+  Q-^     2  cosh  3an/3  =  Q"+  #-*! 
2sinhan0  =  Q'-#/-1J  '   2sinh  3an£  =  #3-#-3J  '  e  - 

where  <2*  =  tan  (£  TT  +  £  i/r,,). 

Thus  (163)  becomes 


m  *,  = 


-1  sin  *  -  <          +     -    sn 

......  (167). 


Writing  u  =  K'm  (163),  *'  =  !»,  and  (165)  becomes 


......  (168). 

(ii)    Use  of  hyperbolic  recurrence  formulae. 

"Write  sin  &  =  sn  08,  *')  =  tanh  ^0',  so  that  (159)  gives 

n  =  -  P  sn2  (iK'  -  1/8)  =  cot*  ^  =  1/sinh2  00',     0<p<£'  ...(169). 
Commencing  with  <£„'  in  (169),  use  the  recurrence  formula 

tanh  (</>'»+!  -  </>,')  =  (6,,/a,)  tanh  <£„'  ............  (1  70). 

We  then  have 

cgS  =  6,72»  =  *',  sn  (y8,  *')  =  tanh  ^»0', 

en  08,  k'}  =  sech  «/>0',       dn  (ft,  k)  =  sec 

oo 

and  -  iZ  (i(3,  fy  =  2  c%  sinh  ^ 

i 
Equation  (160)  then  gives 


As  before,  M  is  calculated  from  the  recurrence  formula 

tan  (<£n+1  -  <£„)  =  (6w/an)  tan  <f>n, 

giving  anw  =  <£*/2R  =  O.    Equation  (163)  then  gives  for  ^  the  expression 
_,  _  sinh  ^'  cos  ^  —  g2  sinh  3^*'  cos  3^>  +  <?6  sinh  5^'  cos  5$  -  .  .  . 
cosh4>'  sin  $  —  ^  cosh  3$'  sin  3$  +  <f  cosh  53s'  sin  5$  —  ... 

......  (172). 

Writing  M  =  ^in  (163),  V  =  %*,  and  (171)  becomes 
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SECTION    X 

NOTE  ON  THE  CALCULATION  OF  THE  THIRD  ELLIPTIC 
INTEGRAL  IN  TERMS  OF  THE  COMPLEMENTARY  A.G.M. 
SCALE 

When  k  is  very  nearly  unity,  it  may  in  some  cases  be  advantageous 
to  make  use  of  the  more  convergent  complementary  scale  («0'=  1,  b0'=  K). 
The  problem  is  to  express  the  various  functions  appearing  in  the  formula 

-  t      ,     .^  sn(a  k}  \     rrr      7\       11       ®  (a  -«, 

n3  (n,  k,  <£)  =  u  +  —  ?  —  j\'   ;  —  7xU^(s  &)  +  |log7T7— 

en  (a,  k)  an  (a,  k)\_  &  ©  (a  +  u, 

.........  (174) 

in  terms  of  recurrence  formulae  involving  the  complementary  A.G.M. 
scale.  The  four  cases  discussed  in  Section  IX,  depending  on  the  values 
of  n,  require  separate  treatment.  The  reader  will  have  no  difficulty  in 
adapting  the  formulae  of  Section  IX  to  the  complementary  scale  should 
the  occasion  arise.  We  may  note,  in  passing,  the  formula 

k'}  .      . 

'" 


which  will  be  found  to  be  useful. 


/a 


IV 


SUMMARY  OF  FORMULAE 

NUMERICAL  CALCULATION  OF  ELLIPTIC  FUNCTIONS 
AND  INTEGRALS 

[Sections  containing  new  formulae  are  denoted  by  an  asterisk  *.] 
First  Elliptic  Integral  «  =  F(+,k)=£ 


Second  Elliptic  Integral   E(<f>,  k)=  S  7(1  -  /P  sin2  <t>)d<f>,  E= 

Third  Elliptic  Integral.  H3  (TO,  k,  <f>)  =  /*- 

/o  I 
Jacobian  Functions. 


sin  <f>  -  sn  ( « ,  k),      cos  «^»  =  en  (u,  k\      *J(l-J?  sin2  <£)  =  dn  (w,  ^t). 
of  Aritkmetico- Geometrical  Means. 


1.   Landen  s  Scale  of  Increasing  Amplitudes 
tan  (<k+1  -  ^>») 

*)=  Lt     1 


I  T7"         ZT'X  /  I7" 
(  ft    —  JCjtiJ\_  = 

^k)-(SjK)F(^tk) 

=  Ci  sin  <^  +  c2  sin  <^  +  ...  +  cn  sin  <£»  +  .... 


32      NUMERICAL  CALCULATION  OF  FUNCTIONS  AND  INTEGRALS 

*2.   Calculation  of  sn  (w,  k),  cn  (u,  k),  dn  (u,  k)  in  terms  of  argument  u. 
Calculate  0n  =  2nanu,  n  being  chosen  so  that  cn  is  negligible.  Obtain 

y» 

0w-i...  02,  <£i>  0o  from  recurrence  formula  sin  (20n_!  -  0n)  =  —  sin0n. 

an 

Then  sn  (M,  £)  =  sin  00,  cn  (w,  &)  =  cos  00,  dn  (u,  k)  =  cos  00  sec  (0!  -  00). 

*3.    jHfo  Hyperbolic  Scale  of  Increasing  Amplitudes 

tanh  (0n+1  -  <£n)  =  (^n/«n)  tanh  ^>n.    Start  with  sinh  <£0  =  tan  <f>. 
)=  Lt  [l/o.. 


—  (G!  sinh  <£i  +  c2  sinh  ^>2  +  ...  +  cnsinh<£n+  ...). 

*4.   Calculation  of  sn  (M,  ^'),  cn  (u,  k'),  dn  (M,  ^')  in  terms  of  argument  u. 

(Iji)  sn  (iu,  k)  =  sn  (u,  k')fcn  (u,  k'),        cn  (iu,  k)  =  1/cn  (u,  k'), 
dn  (m,  A)  =  dn  (u,  k')/cn  (u,  k'). 

Calculate  <£»  =  2nanu,  n  being  chosen  so  that  cn  is  negligible.  Obtain 
<j>n-i  ...  <#>2  >  <#>i  »  0o  from  recurrence  formula 

sinh  (20n_i  -  0n)  =  (cj«n)  sinh  0n. 
Then  sn  (u,  k')  =tanh  00>  cn  (u,  k')  =  sech  00>  dn  (M,  ^')  =  sech  (0!  —  00). 

5.    Calculation  of  K',  E',  and  the  nome  q=e"ttK'lK. 
*  Hyperbolic  recurrence  formula  tanh  Xn+i=  V(a«+i/^n+i)  tanh  2\n- 
F(<t>,k')=Lt  [l/an.Xnn     tanh  Xo  =  N/F  sin  0. 

)Mog^l        .    (Legendre.) 


log  ?  =  5  (£)M-1  log  -^-  -  log  ^  .  (Gauss.) 

1  <*n+l  Ci 

J?"_  1  _  /I    _  l/»2_  .,  2  _  9,,  2  O"-!^  2  _          \ 

—  Jb   —  ^  —  \L  —  -jC0   —  Cj.    —  ^sC2   —  ...  —  a        Cn   —  ...) 


6.   Landen's  Scale  of  Decreasing  Amplitudes 


sin  (2^n+1  -  i^n)  =  -^  sin 
an 

^  — 
an 
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x  —  log  tan  (|*r  +  |</rn)  +  2  to-B,  CB), 

an  2 

2  to    e  >-2g  ten(2^-fr)  |  fi,  tan  (2^-^) 
t  (*"  Cn)  "  2'2  cos  (2^3  -fc)         3  cos  (2fc  -fc)  + 


*7.    Calculation  of  sn  (w,  A;'),  en  (u,  k'),  dn  (#,  A;')  in  terms  of  argument  u. 
Calculate  \j/n  from  sin  $n  =  tanh  (anu),  n  being  such  that  cn  is  negligible. 

Obtain  ^B_a  .  .  .  fa,  fa,  f  0  from  tan  (i/w  -  ^n)  =  ^  tan  <^n.    Then 

«B 

sn  (u,  k')  -  sin  ^0,  cn  (u,  k')  -  cos  i^0,  dn  (M,  ^')  =  cos  (2^  -  ^0). 
*8.   Hyperbolic  Scale  of  Decreasing  Amplitudes 

sinh  (2^B+1  -  «/rn)  =  -5-  sinh  i/r». 
ttw 

^(i/',  >b)  -  l/an  tan-1  (sinh  i^n),        tan  ^  =  sinh  i^0, 

W 

)  ~  -rFW,  *)  =  ^nh  ^0  cosh  (2^  -  ^0)  -  aB  sinh  i^B  -  2'  (^w,  c»). 


2 

*9.    Calculation  of  sn  (u,  k),  cn  (u,  k),  dn  (u,  k). 

Calculate  tyn  from  sinh  \f/n  =  tan  anu,  n  being  such  that  cn  is  negligible. 
Obtain  ^B_j  ...  \f/2,  fa,  i/r0  from  tanh  (^n_!  -  ^B)  =  (cn/an)  tanh  i/^.  Then 
sn  (u,  k)  =  tanh  i/^j  cn  (u,  k)  -  sech  i^0,  dn  (w,  #)  =  cosh  (2^  —  t^0)/cosh  i/f0. 

*10.    Calculation  of  Third  Elliptic  Integral — 4  cases. 

TT  /      7    j\  Sn(a  k)  r     „,      7.       .  ,       ®(a-M,  /l)~l 

ns («,  A-,  <£)  =  M  +  — — =H TT    wZ(a,  A:)  +  i  log  pTT —ft-   . 

cn  (a,  k)  dn  (a,  ^:)  L  ®  (a  +  u,  k)J 

Case  I.   -  k2  <  n  <  0.          Case  II.   —  oo<  AV5  <  - 1.   Hyperbolic  Cases. 
Case  III.    -1<7«<-^2.    (Method  1.    (i)  Circular,    (ii)  Hyperbolic 

recurrence  formulae. 
Method  2.    (i)  Circular,    (ii)  Hyperbolic 


Case  IV.    0  <  TO  <  + 


recurrence  formulae. 


Typical  Solution.   Case  III,  Method  1  (i). 

n  =  -  F/dn2  (a,  k')  =  -  #/(!  -  k'2  sin2  6),     sin  0  =  sn  (a,  k')  0  <  a  <  JT. 

1.  Form  A.G.M.  scale  (a0=  1,  b0  =  k'). 

2.  Derive  an  w  =  </>n/2"  =  *  from  tan  (<£B+1  -  </>n)  =  —  tan  <£„ ,  commencing 

.  ,  a» 

with  <0  =  <. 
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3.  Derive  ana  =  log  tan  (fa  +  frj/n)  from  sin  (2^w+1  -  ^n)  =  —  si 

&H 

commencing  with  ^  =  0. 

00 

4.  Compute  q  from  log  q  =  2  (^)n~1  log  (an/an+1)  - 

5.  Compute  *  from 

sin  2arew  sinh  2ana  +  2g^  sin  4anw  sinh 


_ 

1  +  2q  cos  2anw  cosh  2ana  +  2g^  cos  4aMw  cosh  4ana  +  2g*  +  .  .  .  ' 


(Legendre.) 
6.   Then, 

sin  ^  +  S  (^m,  cn)  -  *], 

2 


sin 


APPENDIX 

In  the  preceding  sections  only  those  formulae  suitable  for  numerical 
computation  have  been  given.  Many  interesting  and  elegant  properties 
of  the  arithmetico-geometrical  scales  in  conjunction  with  recurrence 
formulae  not  specially  adapted  to  numerical,  calculation  are  given  below 
as  problems.  In  some  cases  the  main  steps  of  the  proof  are  briefly  in- 
dicated, and  unless  otherwise  stated  the  final  results  are  new. 

Properties  of  the  A.G.M.  scales  (a0  =  1,  60  =  £')• 

1  .   Prove  that  a»  =  an+1  +  c»+1  ,  bn  =  o»+1  -  c^i  .  Hence  derive  the  series 


....      Gauss,  G.  W.  Bd.  m.  p.  376. 
2.   Prove  that 


and  hence    Vo*  =  *fa0  —  V^  —  Vc4  —  .  .  .  , 

tjon  =  V^+  •Jet-  Vc4-  ...  .    Gauss,  G.  W.  Bd.  ni.  p.  376. 

3.  Prove  that  c2,  +  2c2»+1  =  2  (an  +  o»+1)  cm+1,  and  hence  derive  the 
expansion 

(K-E)IK=  (a0  +  01)  cx  +  4  (a,  +  03)  Cj 

+  16  (a4  +  «5)  cs  +  .  .  .  +  2s*  (a*  +  a^O  c^+1  +  .  .  .  . 

4.  Prove  that  26B_1cn  =  ^e2,,-!  -  2C2,,  and  hence  derive  the  series 


5.  From  the  results  of  Ex.  4,  prove  that 


and  hence  the  identity 

1=«»2  +  22  (2-1-  l)6s_lCj,  +  i  i2»cB2. 

2  0 

6.  From  the  recurrence  formula  of  Ex.  4  prove  that 

2  (6oCl  +  hci  +...)  =  K2  -  1  (ci8  +  c»2  +  c,2  +  -..)> 
and  hence  a.2  =  602  +  ^c03  -  f  (^  +  Ca2  +  c,2  +...), 

«W*  =  «»*-K>-f(^f'<-tf  +  «%*  +  ..)L 

Gauss,  (?.  H7:  Bd.  m.  p.  376. 

7.  Derive  the  recurrence  formula 

log  («»_!/««)  =  log  (bn/bn-J  -  2  log  (b^bn) 

3—2 
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and  hence  prove  that 

\-rrK'  IK  =  log  (4«0/C0)  -  log  (aja^)  -  i  log  (fld/0,)  -  ... 

-(l/2n 
-=  log  (4ft,/e.)  +  1  log  (6  A)  +  !  log  (  W  +  .  •  • 


Gauss,  £.  IF.  Bd.  in.  p.  377. 
[Note  that  \-rrK'  IK  =  (l/2n)  log  (an/cw)  as  n  increases  indefinitely.] 

8.  Derive  the  recurrence  formula 
(1/2"-1)  log  (oH/On-O  =  (3/2B+1)  log  (an/bn) 

+  (l/2"+1)  log  (anlbn)  -  (1/2")  log  (On-A-x). 

Hence  prove  that 

irK'/K=  log  (Itao  W)  +  32  (l/2ra)log  (aw/6n). 

Jacobi,  .P.  N.  §  52,  p.  201. 

9.  Show  that  the  A.G.M.  scale  (»0  =  1,  b0  =  k')  may  be  calculated  from 
the  trigonometrical  recurrence  formula 


commencing  with  sin  00  =  k,  sin  01  =  (1—  #')/(!  +  A;'),  etc. 

Hence  derive  the  following  formulae  : 

(i)  c^^&tani^o,  c2=|-Atan|^0tan^^1,  c3=|-A;tan^^0tan^1tan^^!,, 
Ci2  =  1^  sin  ^i  ,  c22  =  iV^2  sin  0X  sin  ^2}  c32  =  -^^  sin  ^  sin  ^2  sin  Oa, 
from  which,  in  general,  expressions  for  cn  and  cn2  may  be  formed. 

[Note  that  cn/an  -  sin  ^n,  bn/an  =  cos  ^n.] 

(ii)  log  (2JT/*-)  =  Z  log  (1  +  sin  ^)«ST  log  sec  ^0^. 

i  i 


(iii)  log  {*l(2Ktr)}  =  |     log  sec  #„_!. 

(iv)  ^fK'IK=  -  \  log  ?  =  \  log  (4/sin  ^)  -  i  &  log  (1  +  sin  0B+1)- 

[To  derive  (iv),  make  use  of  formula  (61).] 

10.   Show  that  the  A.G.M.  scale  («0  =  1,  b0  =  k')  may  be  calculated  from 
the  hyperbolic  scale  of  decreasing  amplitudes 

tanh  2xn+i  =  tanh2  x»,  or  Xn+i  =  i  log  cosh  2xn, 
commencing  with  tanh  2x0  =  &,  or  xo  =  i  log  {(1  + 

Hence  derive  the  following  formulae  : 

(i)    c^l&tanhxo,     c2  =££  tanh  x0  tanh  Xi,     c3 

C!2  =  ^  F  tanh  2xi  ,  c22  =  TV^  tanh  2\i  tanh  2\2, 

cs2  =  A^2  *ann  2Xi  tawh  2Xs  tanh  2x3, 
with  similar  expressions  for  cn  and  cn2. 

[Note  that  6n/aB  =  sech2xn,  cre/an  =  tanh2x«.] 
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(ii) 


(iv)  \KK'IK=  log  q~$  =  log  (2  coth  Xo)  -  2  (i)*log(l  +  tanhx»«)- 

[To  derive  (iv)  make  use  of  formula  (61).] 

11.   Show  that  the  A.G.M.  scale  («0  =  1,  b0  =  k')  may  be  calculated  from 
the  hyperbolic  scale  of  increasing  amplitudes 

tanh  2i^n  =  tanh2  i/r»-n>  or  «/^  =  J  log  cosh  2i/wi, 
commencing  with  tanh  2^0  =  k',  or  ^0  =  \  log  {(1  +  #)/(!  -  #)}• 

Hence  derive  the  following  formulae  : 

(i)    We.  =  i*"2*".  hJ*»  =  tan11  2<A»,  <W«»=  sech  2^. 


(ii)  log  (2*7ir)  =  £  log  (1  +  sech 


(iii) 

o 

(iv)  \TcK'IK=  log  q~$  =  ^n/2*-1,  as  n  tends  to  infinite  values. 
[To  derive  (iv)  make  use  of  formula  (60).] 

12.  Writing  £,  =  £  {1  +  (1  —  tf-8*")^},  show  that  the  recurrence  formula 
of  Example  1  1  may  be  written 

2^»-n  -  4«A»  =  log  2  +  log  £.  . 
Hence,  denoting  ft  =  (1  -  >/F)/(l  +  \/^'),  derive  the  expansions 

(i) 

(ii)  log  9  = 

(iii) 

(See  also  Example  24.) 

13.  Derive  the  following  series  for  v  from  the  formulae  of  Example  12, 
TT  =  log  x-1  -2a*-  13#8  -  ^x"  -  .  .  .,  where  x  =  |(2*  -  1)/(2*  +  1). 
["Write  k  =  k'  =  I/  V2,  then  Jr=  -ff"  and  log  q~l  =  ».    Note  that 

ar1  =  23-14086,  log  ar1  =  3141602...,  2a^=  '0000070, 
so  that  two  terms  of  the  above  series  give  the  value  of  TT  correctly  to 
the  tenth  decimal.   The  remarkable  approximation  represented  by  the 
first  term  of  the  series  was  first  pointed  out  by  Legendre,  F.  E.  t.  L] 

14.  Given  the  name  q,  compute  the  moduli  k,  k'. 
[The  recurrence  formula  of  Example  1  1  may  be  written 

4f  »  -  2<^+1  =  -  log  2  +  log  (1  -  «- 
=  -  log  2  - 
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Note  that  when  n  is  sufficiently  large  (i)P*l^»»i— JlflgjT*'  We  then 
have 

4iAn  =  2n  log  q-1  -  log  2  -  q*n+\ 

Choose  n  so  that  the  last  term  is  negligible  compared  to  the  first  two. 
The  repeated  application  of  the  above  recurrence  formula  enables  us 
to  calculate  successively  tyn-i  •••  ^2  and  finally  ^0,  in  terms  of  which 
k'  =  tanh  2^0,  k  -  sech  2i/v] 

15.  Derive  the  differential  relations 

11  /  x»\  11  /  ft     \  11  /  \  1  /  \ 

for  all  values  of  n.  Gauss,  G.  W.  Bd.  m.  p.  380. 

16.  From  the  results  of  Example  15  obtain  directly  the  formula 


C 

[In  terms  of  the  A.G.M.  scale  (a0,  b0,  c0)  and  the  complementary  scale 
(«<>'>  ^o'>  co')»  tne  formulae  of  Example  15  give 

and   t- 


But  a0  =  a0',  b0  =  c0',  c0  =  b0'  and  a02  =  b02  +  c02,  from  which  it  is  easily 
proved  that 


Hence  we  conclude  that  \d(^  log  Gf  )  +  —  ^  c?  (—  log  J5-,  )  =  0. 
«w      \z         fw/      «»       \z         fre  / 

We  also  have  the  identity  —,d(-^}+  —^  d  (  —  )  =  0,  from  which  it 
J  an2    W/     an      \an/ 

ft  if 

follows  that  —  log  T^  =  Cn  —,  ,  Cn  being  a  constant  independent  of 

^  On  an 

(an,  bn,  Cn),  («„',  bn',  cn').    In  terms  of  Gauss'  notation  for  the  A.G.M., 
we  note  from  Section  III  that 

«„'  -  M(a0,  c0)  =  2nM(an,  cn)  =  2*0^(1,  cjaj, 
while  an  =  M(a0,  b0)  =  M(an,  bn).    Thus  when  w  becomes  indefinitely 

great,  Cn  =       -  log      =        -  —r-.  log  cf  .  But  if  we  write  k'  =  cn/an, 

,  cn/anj       an 


M  (1,  cn/an)  =  ^Tr/K  =  ^Tr/\og(4lk'~)  as  #'  becomes  small.  Thus,  in  the 
limit  as  n  increases  indefinitely,  Cn  =  -%•*,  so  that 


T      /  1  ,     «TO\     !     an  "I 
Lt   (^S-)  =  ^—" 

W-*.oo  \^  Cn/  «,i      J 


APPENDIX  39 

17.  Making  use  of  the  result  of  Example  16,  prove  that  each  of  the 
ratios  in  Example  15  is  given  by 

V=  Lt   [i7r/(anaB')  .  d  log  («„'/«„)]• 

n-*oo 

18.  Derive  Legendre's  relation  E/K+  E'jK'  -  1  =  ^/(JHT)  directly 
from  the  properties  of  the  A.G.M.  scales  (a0=  1,  b0-k'}  and  («0'-  1,  b0'=k). 

[Differentiating  the  identity  log  (anbn)  -  log  («„_!&„_!)  =  log  («„/£>„)  and 
making  use  of  the  results  of  Example  15,  we  have 

d  log  OA)  -  d  log  (aa_A_i)  =  d  log  (an/bn)  =  2"  (cJIcfldlog  (<z<A)- 
Writing  down  similar  equations  for  «  —  1,  «  —  2,  ...  3,  2,  1  in  place  of 
w  and  adding,  we  find 

d\og(anbn)-d\og(a0b0)^(2cl;i  +  ^+...+2ncH^c0-3dlog(a0/b0)   (i) 
and 


Since  «0  =  aQ'  =  1  and  60'  =  c0,  we  have 

rf  log  (oo'V)  -  rf  log  (aA)  =  d  log  (cA)  =  «o2c0~2  ^  log  (a 

[Ex.  15  with  w  =  0.] 
When  n  becomes  large,  an  =  bH,  so  that 

dlog(an  bn~)-d\og  (a^bn}  =  2dlog  (aB/aw')  =  -  (4/ir)  anan'  c0~*  dlog  (ajbo). 

[Ex.  17.] 

Note  also  the  identity  c0~*d\og  (a0/60)  +  c0'-'2dlog  (a0'/b0')  =  0.  [Ex.  16.] 
Subtracting  (ii)  from  (i)  and  dividing  out  by 

c0-a  d  log  (oo/60)  =  -  c<>' 
we  find 

(4/ir)  anaB'  -  a02 

...  +  2"c»2  +  ...)-  (2C/2 


......  (iii). 

But  we  have  (K-  E]\K=  \  (c02  +  2tf  +...  +  2"cB2  +...),  a«  =  i»/JT,  a0  =  1 
with  similar  expressions  for  the  complementary  functions  in  terms  of 
accented  symbols.  Thus  (iii)  becomes 

*l(KK')  -  1  =  -  2  (1  -  EIK)  +  c02  -  2  (1  -  E'lK")  +  c0/2 
or,  noting  that  c02  +  c0'2  =  1,  we  derive  Legendre's  relation  in  the  familiar 
form 

E/K+  E'jK'  -  1  =  frr/(KK'). 

The  equivalent  of  the  series  (iii)  is  given  by  Gauss  (G.  W.  Bd.  m. 
p.  380),  but  it  does  not  seem  to  have  been  hitherto  pointed  out  that 
Legendre's  relation  follows  immediately  from  it.] 
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19.   Verify  the  following  identities  : 


Gauss,  Nachlass,  G.  W.  Bd.  ni.  p.  447. 

20.   From  the  results  of  Example  19,  prove  that  if  we  form  the 
A.G.M.  scale  commencing  with 


the  wth  term  is  given  by 

an  =  h  (1  +  2#n  +  2#n4  +  ...)2,          6»  =  A  (1  - 


where  log  ^rt  =  2n  log  ^  and  k  =  M(a0,  b0). 

Gauss,  Nachlass,  G.  W.  Bd.  in.  p.  448. 

21.   Prove  that  x  as  defined  in  Example  19  is  identical  with  the 
nome  q  when  («0  =  1>  fy>  =  ^')- 

[When  w  becomes  large, 

=   Lt   si  log  (a.,/0 

7l-*-00   & 

=   Lt   1  log  [(1  -  2^  +  .  .  .)2/4^  (1  +  2*n2  +  .  .  -)2] 

m-*.oo  ^S 

=  -     Lt 


so  that  x  =  q. 

22.   From  the  results  of  Examples  20  and  21,  derive  the  following 
expressions  : 

Jp  =  l 


[From  Example  20  we  obviously  have  a0  =  1,  b0  =  k',  cQ  =  k,  k  = 

23.   If  we  denote  /?  =  (!-  VF)/(1  +  V^')>  derive  the  following  ex- 
pressions : 

!/?  =  (q  +  q9  +  q26  +  ...)/(l  +  2^  +  2^16  +  ...), 

K=  2:r  1  +  2     +  2w  + 
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24.   From  the  preceding  examples  derive  the  well-known  expansion 
=  \p  +  2  ($£)»  +  15  (£#)»  +  150  (&8)13  +  ...  etc.   (See  also  Example  12.) 


q  =  \p  +  2  ($£)» 
25.   If  we  write 


y.  = 
where  log  #„  =  2*  log  g,  prove  that  the  recurrence  formula 

tan  ^,+1  =  (o,+)/^»+1)  tan  2^, 
leads  to  the  result 

Lt 


where  tan  <^0  =  «/#  tan  <£,  and  K  is  given  by  the  ^-series  of  Example  22. 

26.  In  terms  of  Landen's  scale  of  increasing  amplitudes  (Section  IV) 
tan  (</>m+i  —  <£„)  =  (&»/«»)  tan  <£m,  derive  the  following  expansions  : 

dn(u,  £)  =  A(«/>0,  A)  =  a,  +  Ci  cos  <#>!  +  GZ  cos  «^j  +  ...  +c»cos<^)t+  ..., 

00 

^  sn  Menu  =  \  2  2w~1calt_1  sin  «#>»  cos  (2^»,,_1  -  <£„). 

27.  When  (c*+1jan+1)  is  sufficiently  small,  prove  that  <t>*+!  of  Example 
26  may  be  computed  from  the  expansion 

sin  2$. 

i?  sin  4^,  -  i  (c.+^a.^)3  sin  6<^,  -  .... 


28.   Prove  that  the  recurrence  formula  of  Example  26  is  equivalent 
to  An  =  |  (A8,.!  +  IPJI&n-i,  and  hence  show  that 


log  [®  («)/e  (0)]  =  2       log  (O./A.). 

0    « 

[Jacobi,  Fundamenta  Nova  (1881),  §  52,  p.  204.] 

29.   Prove  that  the  recurrence  formula  o»  tan  ^>,»+1  =  A»+1  tan  <f>n  ,  where 
A*  =  ^o*8  cos8  #,»  +  6»2  sin8  <£,,),  may  also  be  written  in  the  forms 

sin  <^,1       2o,+1 


0.  +  A,5      tan^,,  "   Vo^   Va.  +  A, 


30.   Prove  that  as  a  result  of  the  recurrence  formula  of  Example  29, 
<^B+1/Aw+1.   Hence  derive  the  following  formulae  : 

,  k)  =   Lt  (^n/a,),  where  </>  =  <^0,  sin  <^0  =  sn  («,  4), 
»-*•• 

=  2<a  cos  ^>0  sin  <^  +  4c,  cos  </>!  sin  ^>2  +  .  .  .  +  2*cm  cos  ^,-j  sin  ^  +  .  .  .  , 

log  [e  («,  *)/e  (o,  *)]  =    a-  log  [ao^/Co.  +  A,)]. 
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[The  recurrence  formula  which  lies  at  the  basis  of  the  results  obtained 
in  Examples  29  and  30  was  originally  given  by  Gauss  in  1818  in  his 
paper  "  Determinatio  Attractionis,  etc."  G.  W.  Bd.  in.  (1866),  p.  353, 
and  was  further  discussed  by  Jacobi,  "Zur  Numerischen  Berechnung 
der  Elliptischen  Functionen,"  G.  W.  Bd.  I.  (1881),  p.  351.  It  is  easily 
seen  that  these  formulae  are  not  adapted  to  rapid  numerical  calculation.] 

31.   If  sin  <£0  =  sn  (u,  fc)  prove  that 

tan  anu  =  —  *    2  '"  —  ^^-  tan  <£0  . 


[The  result  follows  from  Gauss'  recurrence  formula  of  Example  29  ; 
see  also  Jacobi,  F.  N.,  §36,  G.  W.  Bd.  i.  p.  154.] 

32.  Prove  that  Landen's  recurrence  formula 

tan  fyn  -  </rn+1)  =  (cw+1/a».,.i)  tan  i/^+1 
may  be  written  tan  2x»+i  =  /s/(cw/a»)  tan  Xn  if  we  write 
tan  \l/n  =  VCon/c.)  tan  x». 

33.  If  we  make  use  of  Landen's  scale  of  decreasing  amplitudes 

sin  (2\ffn+i  -  i/O  =  (&„/«„)  sin  $n, 

commencing  with  i/^  =  |TT,  prove  that  K'—  (2K/ir')  log  tan  QTT  +  |^TC)  when 
n  increases  indefinitely.  Hence  show  that  the  nome  q  may  be  computed 
from  the  formula  q  =  cot2  QTT  +  £</^)  =  (1  —  sin  «An)/(l  +  sin  "/0-  Show, 
furthermore,  that  Legendre's  relation  leads  to  the  identity 

an  (1  -  sin  «/O  =  (frr/ST)  {1  -  tanh  (irK'^K)}  =  2  fyn,  cn), 

2 

where  2  ($n,  c»)  is  the  series  of  equation  (75). 

2 

34.  Prove  that  if  the  A.G.M.  scale  (a0,  b0}  is  computed,  the  circum- 
ference I  of  an  ellipse  of  semi-axes  (a0>  bo)  is  given  by 

I  =  (27r/«n)  (a?  -  2c22  -  4c32  -  8c42  -...)• 

35.  If  we  denote  by  r-^  and  r2  the  greatest  and  least  distances  re- 
spectively of  a  point  on  one  circle  to  the  circumference  of  a  coaxial  one, 
show  that  in  terms  of  the  A.G.M.  scale  («0  =  n  ,  b0  -  r2)  the  coefficient  of 
mutual  induction  is  given  by 

M=  {W/an}  [c?  +  2c22  +  4c32  +  .  .  .  +  2ft-1cn2  +...], 

or,  when  the  circles  are  very  close  together,  in  terms  of  the  comple- 
mentary scale  (a0'  =  ri  +  r2,  b0'  =  rl  —  r2),  by  the  formula 
(r12-2c2'2-4c3/2-... 

~~    - 


,r( 

2ira» 


[King,  L.  V.,  Proc.  Boy.  Soc.  A,  Vol.  c.  1921,  p.  63  et  seq.    Several 
numerical  illustrations  are  given  in  this  paper.] 
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